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Abstract 

This article considers, along with its various direct extensions and ramifications, a technique for 
analytic interpolation over the exterior of a unit disk using complex poles in the interior — as well as 
corresponding techniques for the exterior of a real unit disk and for the interior of a real and complex 
unit disk. More specifically, a simple closed-form technique is presented for approximating analytic 
functions in the exterior of a complex unit disk by linear combinations of simple poles located in the 
disk's interior at distinct points. This is accomplished by developing special kernel spaces labeled dual- 
access collocation-kernel spaces. In particular, it is shown that for two common norm settings a bounded 
point replicating kernel exists with the ubiquitous form l/{z — z^), where z is in the exterior region of 
interest and where the fixed source point Zk is inside the disk. When the ZkS are specified, both norm 
settings yield exact closed-form equation sets for the associated pole strengths. The two norm settings 
are the standard inner product one for functions prescribed on the unit circle (where boundary value 
differences are minimized) and the Dirichlet integral one (where field energy differences are minimized 
over the whole exterior of the unit disk). Besides 1/(2: — z^), logarithmic point source kernels are also 
considered in both the real and complex planes. Higher order poles are also considered. Results of 
numerical tests are given. Certain educational possibilities are addressed. Relationships to the Cauchy 
integral formula, Bergman kernel theory and Szego kernel theory are also addressed. For one of the 
norm settings — Dirichlet-integral dual-access collocation space (DIDACKS) — this is the first of several 
planed articles and these other articles will address various generalizations and applications of the basic 
technique presented here. 
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1 Introduction 

First consider the paper's content and structure. This article develops a simple closed-form interpolation 
approach that can be used to approximate complex analytic functions in the exterior of a unit disk by a finite 
set of spatially distributed poles in the disk's interior. This approach rephcates values of the analytic function 
to be fit (and/or its derivatives) exactly at specified points and thus requires only limited information at 
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a finite number of selected exterior points. Furthermore, this closed-form point replication approach also 
minimizes the difference between the given function to be approximated and the approximating form (i.e., 
the modeling difference) for one of two common norms of interest and it produces a set of exact closed- 
form linear equations that can easily be solved for the corresponding Nk pole or point source strengths. 
The dual interpolation and approximation nature of the solutions obtained from these equation sets means 
that they are actually interpolation fits. Mathematically this is accomplished by developing special kernel 
spaces that, while closely related to reproducing kernel Hilbert spaces, are different from them. (Notice that 
since it is assumed that noting is known about the function to be modeled in the interior of a unit disk, 
the strength of the individual poles cannot be found in the usual way by performing interior line integrals 
and then applying the theory of residues.) The general kernel structures associated with these spaces are 
briefly outlined in Section [5] and various associated notational and nomenclatural conventions are given in 
Section [31 This is followed up with a discussion of the approximation space setting in Section 2) Although 
the corresponding theory for functions is not explicitly considered until Section [51 in order to properly 
set the stage for their introduction the foundation for their use is also laid in Section [JJ hence, there are 
two primary classes of functions that must be dealt with here: (a) A restricted class of C functions that are 
analytic over the exterior of a complex unit disk, (b) A restricted class of functions that are harmonic 
over the exterior of a real unit disk. These two classes of functions are defined in such a way as to be 
coextensive (i.e., in direct correspondence). In the end, this is done by requiring that the class of admissible 
analytic functions be restricted to those that have a power series representation of the form ^ and then 
by assuming that corresponding functions can, in principle, be obtained through a process of harmonic 
completion. While other ways of defining the relevant classes of admissible functions exist, one reason that 
this approach is used is that it reduces the level of the core material (basically the material from Section[5]on) 
to a bare minimum. The accessibility of this core material, in conjunction with certain other factors, opens 
up several unique educational possibilities that are partially outlined below and in Appendix A. Section [2| 
starts out by considering the class (b) functions just mentioned, but defined over a more general domain, 
which corresponds to a real functional analysis setting with the functions defined over some sub-domain (and 
possibly an unbounded exterior region) of R'^. In Section [2] it is shown that the direct analytic/real function 
correspondence mentioned above side-steps standard Hilbert space and Sobolev space theory in favor of a 
structured (i.e., augmented) pre-Hilbert space setting. As partially discussed in Section[4l this has the effect 
of circumventing certain issues related to the kernel forms used, as well as other issues tied to matters of 
physical interpretation. Section [4] also addresses various other foundational topics such as the uniqueness 
of simple pole expansions (and real dipole expansions) and the possibility of interpolation fit solutions that 
have good modeling properties while being the solution of an ill-conditioned equation set. At the end of 
Section [4| a transition is made to the more concrete mathematical matters that frame the core of the article. 

Section [5| gives a straightforward and self-contained example of these kernel structures within the context 
of the first norm setting. This first norm minimizes the modeling differences on the unit circle and is formally 
proportional to the standard Hilbert space norm of Lebesgue square integrable functions over this circle. 
This norm is labeled the standard norm so as to distinguish it from this regular norm and its attending 
Hilbert space backdrop. The space associated with this first norm is called the standard-integral dual-access 
collocation-kernel space (SIDACKS). (Since acronyms are common throughout the paper a list of them is 
included in Section [31 for the reader's convenience.) Connections to the Cauchy integral formula, Bergman 
kernel theory and Szego kernel theory are also addressed in Section [5l The second norm is introduced in 
Section [6] and it is based on the well known Dirichlet integral, which means that it is an energy norm for the 
modeling difference. The associated space for this inner product structure is called a Dirichlet integral dual- 
access collocation space (DIDACKS) and the primary mathematical relationships for it are also presented 
in Section [HI 

Up to Section [7| the article is primarily concerned with analytic functions in and on the exterior of a 
complex unit disk (i.e., f{z) for z S C and \z\ > 1). As noted above, the general perspective adopted 
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here is that problems in C and can largely be considered interchangeable since: (a) a function that 
is harmonic over some region of interest can be uniquely completed to form an analytic function / 
over the same complex region (c.f., Definition [71 1 and surrounding discussion) and (b) both the real and 
complex components of an analytic function / can separately be considered harmonic functions (where 
the complex component is usually simply ignored). Thus the technique presented in the paper also forms 
a basis for framing associated harmonic fitting problems in where logarithmic sources commonly occur. 
Section [7] first extends the technique to handle complex logarithmic point sources and then these complex 
plane results are carried over to the real plane, M^; hence, closed- form equation sets are shown to exist for the 
strength of logarithmic (or dipole) point sources in both norm settings. (In Section [7] it is pointed out that 
it is actually more efficient to map the norm structures themselves from C to K.^ rather than use the above 
harmonic to analytic mapping and then the return mapping.) Next extensive, but not exhaustive, numerical 
testing is performed in Section[5] Given the point replicating nature of the algorithms here, these tests serve 
not only to check the software implementation, but also, to a certain extent, the formalism itself. As an aside 
it is worth noting that educators looking for interesting examples of ill-conditioned matrices should consider 
the results presented in Section [8l but by the same token those who are inexperienced should exercise a great 
deal of caution when dealing with systems that have condition numbers as large as those that are considered 
in Section [8] — expecially if any computational or measurement errors are present. 

While the focus on the exterior of a unit disk is retained up through Section[8l in the next and final section 
(Section [HI) it is shown that the basic approach can also be developed for interior approximation problems, 
where the region of interest is the interior of a unit disk and the approximating poles are located outside the 
unit disk. This means that all of the various embellishments to the formalism that were developed in the 
prior sections can also be carried over to this interior setting. Appendix A contains additional discussions of 
certain side issues: In addition to the educational factors mentioned above. Appendix A also considers certain 
pedagogical aspects of the material in the main body of the article. Finally, Appendix B briefly addresses 
the side issue in the complex setting of what specific expression to consider for the Dirichlet integral when 
the functions are either analytic or antianalytic. Before proceeding readers may wish to note that kernels are 
generally placed on the left and that the complex conjugate of the first factor is taken in an inner product 
rather than the second — in Section [3] the acronym left conjugate form (LCF) is introduced to define this 
convention. 

Next several issues pertaining to the paper's notation and style will be addressed, but first, however, it 
is useful to briefly consider the paper's origins. The method described herein had its genesis in a relatively 
recent detailed reexamination of gravity modeling linear least squares (LLSQ) and non-linear least squares 
(NLLSQ) techniques that the author first developed and successfully applied about a quarter of a century 
ago at the Naval Surface Weapons Center Dahlgren Division and, as such, the method should be considered 
technology spin-off. As noted in the last line of the abstract, the associated general interpolation/modeling 
framework, as well as certain applications settings and theoretical implications, will be described in a sequence 
of articles, of which this is the first. Most of the articles in this sequence are already largely complete. Even 
though they are labeled as a sequence, most of these articles should be more-or-less independent of the 
others and so, except for a few exceptions, can be read out of sequence. Later in this sequence a brief 
discussion of the general goals of this sequence is planned. After the core articles in this cycle are complete 
an overall survey article will be written that may contain this discussion of goals, provided that time and 
circumstances permit. Finally the author wishes to thank his group leader, Larry D. Wilkerson, and his 
branch head, Howell J. Boyles, for their patience and forbearance during the gestation of this sequence 
of articles and, in advance, for their review or perusal of each of the articles in the sequence. As usual, 
any lapses of judgment with regards to content and style, as well as the occurrence of typos, are the sole 
responsibility of the author. Due to the overall accessibility of these articles, since mathematical errors or 
typos can be easily detected from the surrounding context and since DIDACKS implemented algorithms 
display an underlying point collocation (or interpolation) property, which serves as a self-consistency check. 
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typos are not considered to be a significant issue in themselves, but by the same token a reasonable effort 
will be undertaken to eliminate them. 

Having considered it structure, content and origins it is useful to examine the article's overall style. The 
detailed discussions that ensue from this examination will occupy the rest of this section. Due to educational 
and other factors, in part of the discussions that follow certain topics that are not normally addressed in 
research articles are addressed and this results in some encroachments of the nominal research article style 
conventions. To otherwise limit the extent of these stylistic encroachments, these discussions occur here 
primarily only in the rest of this section and in the appendices of this and future articles. 

Given these, as well as other style and notational issues, the reader may wonder why the paper was written 
as it is. First, while it may not be clear to the reader from the content of this article alone, a somewhat 
irregularly shaped hole was perceived by the author to exist in the literature that entailed pure, applied and 
educational factors and, after attempting to balance all of these aspects, the approach initiated here seemed 
to be the most appropriate way to undertake an initial effort at partially filling this hole. Second, either the 
material here is sufficient to stand on its own merits or it is not. If it is, then it can endure these stylistic 
abuses. If it cannot, it will be ignored in any case. Aside from the overview article dealing with goals, the 
general stylistic tack taken in the other articles in the sequence will be that of a nominal research journal 
article where it is entirely up to each individual reader to discern whether the material in each article does or 
does not have intrinsic merit with regards to his or her own framework of research or educational interests. 
While this nominal journal style could have greatly simplified matters here, it would have been necessary 
to abandon all educational factors (in the end, these factors often do not take care of themselves so they 
must be taken into account at the onset). Those readers who are not otherwise interested in these larger 
issues — namely educational aspects, aspects dealing with the overall setting of the approach and an attempt 
at assessing the novelty of the technique — may simply skip the rest of this section (as well as Appendix A) 
and proceed directly to Section O 

With regard to the issue of novelty or uniqueness of the approach, which is discussed at length at the end 
of this section, in the end it seems premature to make a definitive assessment here and it is the hope that the 
passage of time will allow for a better perspective. At present what seems to be mainly lacking, aside from 
those connections pointed out in the text, are known primary references — if any such references exist at all. 
Note that while there may be some truth in the unthinking reflex reaction of some that anything of interest 
in basic elementary complex variables that can be done has already been done, this applies only to the extent 
that this well known material is looked at from a more-or-less standard perspective; however, if this standard 
material is looked at from a sufficiently different point of view new results may indeed be possible. What 
seems to be different here is to take the idea of a dual-access collocation kernel (DACK) and its ramifications 
seriously. At the very least, this concept seems to shed some new light on otherwise well known results. It 
also introduces a coherence and unity into otherwise disconnected results in complex variables and potential 
theory. Finally, in some real sense it is this overall context that transforms a collection of interesting but 
isolated results into a potentially useful tool kit. Setting forth this context is one of the overall goals of the 
sequence of articles. 

For better or worse this article can be viewed from either of four different perspectives: 

{ I } As a research article. (Or more specifically one on analytic functions and the theory of poles and 
related topics in M^.) 

{ II } As an article addressing educational aspects. 

{Ill} As a review article. 

{IV} As the first in a sequence of articles that introduce either an author's unique line of research or an 
author's unique take on existing material. 
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It is not unusual for various articles to simultaneously display three of these attributes. In fact, most research 
articles contain a brief introductory section that sets the new results in context by introducing and briefly 
surveying the relevant research literature, thus encompassing attributes {1} and {III}. (As noted below, 
the literature is not surveyed here since, aside from long established results that are discussed at length 
later, no closely related research results are at present known by the author to exist.) By the same token, 
a research article that is the first in a sequence may quite often simultaneously encompass some mixture of 
attributes {I}, {III} and {IV}. Likewise, the first of a sequence of educational articles might be expected 
to encompass aspects {II}, {HI} and {IV}. What is unusual here is that, to some degree, all four aspects 
are simultaneously displayed. This is unusual since attributes {1} and {II} are generally considered to be 
diametrically opposed. Thus, the usual strategy that is employed by most researchers is to simply trump 
any significant educational (and review) possibilities and produce a pure research article. As discussed below 
and in Appendix A, this was not done here since these educational prospects are thought to be significant 
and unique here and the overall impact on the research aspects of the article that were incurred by taking 
these prospects into consideration was thought to be acceptable (but, of course, in the end this is a matter 
of personal taste and the same trade-offs will not necessarily be made even in the future articles of this 
sequence). This does, however, lead to certain inevitable complications of style and notation and it is thus 
useful to consider these complications from each of the four perspectives just listed. 

{1} This article is primarily intended as a research article, but when viewed from the perspective of 
any one of these four article classes, there are several noteworthy complications here due to the concurrent 
educational and review aspects of the article. First, the article is longer than it might otherwise have been, 
but in the main body of the article the extra steps that are included in some derivations, as well as the extra 
explanations, should lessen rather than increase the effort required by researchers to read and understand 
it. Otherwise, when the results are considered on their own terms, a fairly mathematically direct approach 
results. Second, from the perspective of a research article on analytic functions and the theory of poles, 
there are many various ramifications and implications that are not brought up, so upon closer examination 
the article has an unsatisfactory open-ended nature to it. (A few of these possibilities are addressed below, 
but most are simply by-passed.) Third, also from the perspective of a research article on analytic function 
theory, some non-standard notation is used. In particular, as previously noted one convention adopted 
here is that kernels are usually placed on the left and the complex conjugate of the first factor is taken 
in an inner product rather than the second factor. The reasons for doing this are discussed at length in 
Section [3l where it is argued that here there are strong pedagogical advantages inherent in adopting this 
convention, which include directly leveraging analogies from linear algebra and operator theory, as is done, 
for example, in the theory of integral equations. Introducing kernel based approaches quickly and painlessly 
is very important, not only for areas of pure mathematics, but for various diverse areas of applied science. 
Aside from providing a minor nuisance to researchers, there is little real down-side here to adopting this 
inner-product conjugation convention, provided that the use of this convention is largely confined to an area 
that can be considered a new specialty or distinct sub-discipline, which is assumed to be the case here. It is 
very important that conventions not be needlessly mixed in long established areas of mathematical research, 
since added complexity and confusion may otherwise result. 

{II} The educational possibilities just alluded to are worth further elaboration, if even from a limited 
preliminary perspective since, as just noted, they have a bearing on the style and notation of the paper 
itself. In fact, as discussed at the end of Appendix A, the acronym DIDACK was chosen not only to imply 
direct connections to Dirichlet's integral, but to imply certain educational and self-educationally attributes 
of the basic approach — especially in the self-educational sense of an autodidact. The other main educational 
possibilities are also briefly listed in Appendix A, but any actual realizations in this area will have to depend 
upon the efforts of others since the author plans to make no systematic efforts in this direction. With 
regard to educational factors then, the goal of this series of articles is merely to preserve the possibility of 
educational aspects in case others wish to follow up on them. Moreover, besides notation, the only systematic 
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pedagogical concessions that will be made in future articles is that some effort will be made to keep them 
reasonably accessible. 

Finally, there is a problem from an educational perspective with the article that is common to the three 
other perspectives also. Specifically, in a graduate class room setting, it is incumbent on lecturers to point 
out what is standard known material and what is novel and unique — unique both with regards to what has 
been recently done by others and what the lecturer has added to it. This is necessary since students are 
often at a loss to ascertain the boundary separating these different types of materials for themselves and one 
expectation of a core course in complex variables is to provide practitioners not only with an understanding of 
standard core concepts, but with a common shared language so that they can communicate with others in a 
professional capacity and understanding what others know or do not know is a key part of this communication 
process. The problem here is that even researchers who are very experienced with the standard core material 
of complex variables may well vary on their assessment of what is actually new and unique here. It is for 
this and other reasons that a fairly long discussion of this topic is included at the end of this section. 

{Ill} Next, consider problems with the article when it is viewed from the perspective of a review article. 
First, and primarily, a review article of a mathematical field must be about some established area of research 
or sub-discipline and it generally should contain extensive references to the known work. DIDACKS theory 
is not an established research area yet so there are no appropriate references; moreover, at present there is 
very much an open-ended quality to the research presented here. Nevertheless, the material presented here 
forms a novel perspective and sunniiary of much of the standard fare of complex variables, but if the reader's 
intent is simply to review complex variables, then there are many fine existing texts. 

{IV} Since this article was expressly written as the first in a sequence of research articles, it perhaps 
works best in this role, but there are still issues. First, the overall goals are not addressed here. Second, 
no attempt will be made here to indicate applications, but applications related issues will be addressed in 
future articles of this sequence. Third, with regards to pure mathematics research, although absolutely no 
attempt will be made here to convince the reader that other actual additional research possibilities exist, 
there are three general classes of research areas where DIDACKS theory may, in some sense or other, be 
suggestive: 

Three General Types of Pure Mathematics Research Areas Worth Considering 

(A) Well established and otherwise "mined-out" areas, which are worth re-examining from the perspective 
of the approach presented herein. 

(B) Currently active areas of research that may yield otherwise unexpected results when examined from 
the perspective of the approach presented here, in the sense that DIDACKS theory may prove to be 
suggestive to an otherwise receptive researcher. 

(C) DIDACKS theory may, in some way or other, suggest entirely new lines of research. 

To make all of this slightly less abstract, an example from each of these research areas that might be worth 
considering from a DIDACKS point of view is given next: 

An Example of Each Type of Research Area and DIDACKS 

(a) There are classic theorems in the theory of analytic functions and poles that link the value of an 
analytic function to its value at a sequence of points. There are also theorems in approximation theory 
that give error bounds and that link the value of a function to a series of partial sums of approximating 
functions. DIDACKS theory deals with both of these classes of results and thus may imply that there 
are some unrecognized links between them. 
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(b) When considered from the DIDACKS perspective, there may be interesting insights that emerge from 
analytic matrix polynomial theory. (For classic references to analytic matrix polynomial theory see 
those in [3] and [TTj. [T7] also gives a self-contained overview.) 

(c) The concept of dual-access kernels and dual-access collocation kernels are introduced here and they 
may be worth exploring in their own right. 

These particular examples are given primarily in the hopes that others will, at least, consider them in 
passing. First, present plans are to not explicitly address either (a) or (b) further here or in any other articles 
in the sequel. Also with regards to type (B) theories, there are many other approaches in analytic function 
theory that apparently have little in common with the approach taken here, such as Carleson interpolation 
and Pick interpolation (or Nevanlinna-Pick interpolation), but that are always worth examining from a 
different perspective due to their inherent richness and one way to do this would be to compare and contrast 
them with the DIDACKS approach outlined here. (For example, with regards to Nevanlinna-Pick theory, 
in spite of overall differences, there are certain obvious research questions that stick out here and are worth 
exploring.) Second, item (c) frames part of the discussion of kernel structures given in Section[2]and may be 
partially addressed in future articles in the sequence; nevertheless, even item (c) will not be, in any sense, 
be completely settled in this sequence of articles. 

Finally, the motivation behind item (c) is worth briefly considering. To date, only symmetric reproducing 
kernels have received the full attention of the mathematical community. The conventional wisdom seems to 
be that since an asymmetric kernel can generally be reformulated as a symmetric one and symmetric kernels 
have very nice properties only symmetric kernels are worth considering; however, specific counterexamples 
contained in this article indicate that this may not be a universally valid assessment. The hope is that the 
discussion in Section [51 when coupled with other preliminary analysis and examples given in later articles in 
this sequence, will convince some researcher (s) that dual-access kernels and dual- access collocation kernels 
are worth exploring in their own right. The approach to this issue taken here is to try to place SIDACKS 
and DIDACKS theory within a larger nomenclatural context, which requires drawing certain distinctions 
that are not usually drawn — these distinctions are taken up in Section[2l The nomenclatural complexity that 
ensues partially results from the fact that the perspective taken here is somewhat unique and partially from 
unavoidable conflicts with standard kernel naming conventions; however, as indicated above, these linguistic 
complications are offset by an easy overall accessibility of the article itself so that relevant implications can 
be readily drawn from the surrounding context. 

What is new here? 

Generally this is a question asked by readers who view mathematical progress as an exercise in taxidermy 
rather than as a growing dynamic enterprise and such readers are not usually interested in the possibility of a 
new synthesis or framework; nevertheless, this question is worth exploring in some detail at some point in this 
sequence of articles so this analysis is undertaken here. Although this same question can be entertained again 
at various points in the future and for other DIDACKS contexts, the general points raised here should also be 
relevant in other settings as well so this type of analysis will not be repeated. [Likewise, while the emphasis 
in this analysis is on the overall kernel setting, one can raise the same question with regards to various other 
specific (and obvious) concepts (such as the paired point formalism), which in some sense or another all 
appear to be new, but such items will not be discussed separately.] The convention of research articles is 
followed here and in other articles in this sequence that assumes all results, when not otherwise noted, are 
new unless they are obviously known results. In this connection, a concerted effort was made to search the 
relevant literature to ascertain significant precedents, but aside from several obvious connections to Szego and 
Bergman kernel theory that are discussed later in the text, there appeared to be none. As discussed below. 
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this is very surprising in itself, but it would be even more perplexing if the general approach developed in 
the sequel had occured somewhere or other in the literature and had been generally ignored, since this would 
otherwise seem to entail a significant oversight on the part of the broader mathematical community itself. 
Although suitable references are not known, it does seem very likely that in the (mainstream?) literature the 
fact that Szego and Bergman kernels can be recast in the form l/(z — z^) has long been recognized; however, 
without the concept of a dual-access collocation kernel and the recognition of its general significance this 
fact was dismissed as a mere curiosity. It would seem, then, that what is obviously new here is the overall 
perspective and with regards to particular results there is some question as to where the actual boundary lies 
between what is known and what is new. However, even this assessment is not entirely clear, since there are 
at least six inherent difficulties in actually performing a complete and exhaustive search of all the relevant 
current and historical literature to ascertain precisely what is actually known and what is not; moreover, it 
is not even clear whether this assessment is too conservative or not conservative enough. 

First, there is the sheer magnitude of the current mathematical literature that would have to be searched. 
For example, it might initially seem that only the complex setting (including analytic function theory and 
its various offshoots) would have to be searched since results that are known for real disk geometries are 
invariably carried over to the complex setting, but even this assumption need not be true since the work here 
is also tied to other broad areas of mathematics, including approximation theory, potential theory, integral 
equation theory and functional analysis (specifically reproducing kernel Hilbert space theory), which each 
have their own independent literature. Second, a broad search of the current literature on analytic function 
theory and complex variables is more than sufficiently daunting, but matters are somewhat worse than usual 
here since most of the topics discussed realistically could, and perhaps should (and thus may), have been 
developed 50 years ago or much more. This means that the literature to be searched must potentially span 
the entire historical record of a better part of a century. Third, this search is compounded by the fact 
that various countries (including all of the major historically industrialized countries) have their own rich 
and varied traditions in complex variables and that the literature of each of these countries has their own 
associated technical and national languages. Fourth, there are many separate disciplines within each of 
these countries that practice complex variables (i.e., applied mathematics, physics and various engineering 
subspecialties). Fifth, the difficulty in assessing what is new here is also greatly compounded by the fact 
that many of the fundamental relationships obtained in the sequel are derived by well known elementary 
techniques. [For example, the derivation of (|14p uses several ubiquitous series representations along with 
the correspondence z ^/z; moreover, itself is an alternative form of the Cauchy integral formula.] 
Clearly, when a result is obtained by elementary means it is tempting to assume that it is well known when, 
in fact, it may not be, but at any rate it clearly has the unfavorable by-product of greatly extending the 
range of literature that must be searched to include literature with a more elementary scope. Sixth, on closer 
inspection the material presented here is somewhat broader than one might expect and this induces distinct 
literature assessments issues of its own. For example, even for a well trained historian of mathematics there 
is usually a subjective element involved in determining all relevant precursors and their impacts on some 
broad topic in mathematics; moreover, it is generally recognized that the task of assessing origins is especially 
difficult and often fraught with potential controversy. In summary, all that can be asserted in the end is 
that, besides Szego and Bergman kernel theory, no direct precedents were uncovered by the author in his 
literature searches. 

Temporarily setting aside connections to Szego and Bergman kernel theory, the lack of obvious direct 
precedents in itself leads to a quandary. First, given that the results presented here seem to have a fundamen- 
tal bearing on the theory of poles, if something along very similar lines had been obtained sometime in the 
past and its significance had been overlooked by the broader pure and applied community, this in itself would 
have been very surprising — in fact, much more surprising than the fact that it had not been discovered at all. 
A matter of judgment, however, does enter here though: if the reader perceives the content of the present 
article as being very specialized and of limited interest then since many disconnected specialized results exist 
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in the literature it might be reasonable to assume that some significant precursors exist in the literature that 
are unknown to the author. (In this regard, it is not unusual for specialized results to be perceived by a 
researcher as new when they were, in fact, originally done 50 or more years ago.) But even here there is a 
problem since the judgment that some results are of limited interest must be almost universally held or else 
those who have a different opinion will integrate them into their own research efforts or expository writings. 

In summary, it seems likely that there is a real quandary here since otherwise how could something so 
simple have been overlooked, but there is a simple resolution of this paradox. To the extent that lessons 
from the history of science also apply to pure and applied mathematics — and there are a sufficient number of 
concrete examples from the history of mathematics to argue that this is indeed the case — when a new point- 
of-view emerges different questions are asked and entirely new results emerge. It is the author's contention 
that the DIDACKS programme can be considered a fresh point-of-view, but the breadth and full extent of 
this point-of-view may not in any sense be apparent until this series of articles is complete. 

Finally, it is also useful to look at the content of the present article with the express goal of placing 
some overall bounds on what the original content actually is here. Thus to set a crude upper bound on 
the originality of the work here first consider the concrete links between the formalism presented here and 
Szego and Bergman kernel theory, which are discussed at length in later sections. These links lead one to 
think that specific replication properties [such as (fT4|) and (j44p ] should probably not be considered new, 
which immediately gives an upper bound on the possible originality of much of the work presented here. 
Alternatively, observe that for C and disk geometries Szego's and Bergman's kernel are, in some sense, 
equivalent to Cauchy's integral formula [as are all similar integral results such as Poisson's integral or Green's 
theorem] and they are thus themselves derivative results when considered at least from a sufficiently basic 
elementary level; nevertheless, each of these kernels is of great separate interest in itself, not only for its basic 
properties, but for the copious related research that it has spawned. In this sense, one might argue that the 
true measure of originality and significance can only be made in historical hindsight. [One might argue that 
the recast form of these kernels given by and have some small claim to distinctness thus calling 
this bound into question.] For the lower bound consider the following. In the potential theory literature 
there seems to be no natural way of performing simple logarithmic or dipole fits and from this void alone one 
can argue that the approach initiated here helps to fill a theoretical hole. Moreover, there seems to be no 
good elementary formalism for performing simple pole fits — which can be considered a natural extension of 
elementary complex variable theory. From these deficiencies one can (tentatively) conclude that at least the 
overall approach presented here is in some sense new, or that at least within the context of these applications 
it is new, thus giving the desired bound in the other direction. 



2 Basic Kernel Setting 

In order to frame the basic general kernel setting of the DIDACKS approach it is useful to compare and 
contrast it with standard reproducing kernel Hilbert space (RKHS) theory. RKHS theory has numerous well 
known merits p^, HP, ^25] . It is useful to introduce RKHS conventions by way of a conceptualized example. 
For real valued functions defined over the real interval X := {x G [a, b] \ a, b and x € M with a > 5 } 
and with the inner product (/, g) :— Jj- f{x) g{x) dx, a reproducing kernel Kr has the property that 
{Kji{x, y), f{x)) = f{y), where x, y 2 and Kn is assumed to be bounded and it is symmetric: Kji(x, y) = 
Kji{y, x). Reproducing kernels for other inner products and other subregions of C and M" are defined 
analogously for n >2. These subregions fl are simply the allowed domains for the functions of interest. For 
a fixed subregion and associated inner product it is easy to show that when a reproducing kernel exists it is 
unique. This uniqueness underlies the usual naming convention for RKHSs, which generally attach a specific 
mathematician's name to an associated reproducing kernel. A compact and precise alternative description 
of reproducing kernels can be based on the Riesz representation theorem ([H], pp. 242-244 of Ch. 9). 
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Two RKHS examples of particular import follow immediately as generalizations of the above example. 
First, the class of reproducing kernels where the norm is defined as a straightforward integral of |/p over fi 
are known as Bergman kernels. Secondly, when the unit disk interior in C is the region of interest and the 
norm is the integral over the unit circle of the square of the magnitude of some complex valued function, 
then the associated kernel is called the Szego kernel. As noted elsewhere the Bergman kernels for the interior 
and exterior of a unit disk and the Szego kernel will prove to be of special interest since they are closely 
related to the kernels studied here. 

When Bergman and Szego kernels are put in a general context one arrives at the left hand side of the 
following nomcnclatural hierarchy structure: 



which proceeds from general to particular on both the left and right hand sides. Spaces of roughly the same 
generality are displayed on either side of each line. The immediate goal of the remainder of this section is 
thus to explain the right hand side of this nomcnclatural structure. 

Before introducing the concept of a dual-access kernel (DAK) it is expedient to clarify the nature of 
allowable subregions. While the applied parts of the paper from Section [5] on center on disk geometries for 
where il C C or C K^, in this section not only are more general geometries in C and considered, 
but also often the cases ft C M". (Here and in the sequel for finite n > 2 will be denoted simply by 
M".) For subsequent developments it is expedient to impose additional structure on Q. In particular, it is 
assumed that C or M" is partitioned into two regions: il and its compliment fi', where at least one of these 
two regions is simply connected. Further it is assumed that O contains the common boundary separating 
these two regions (denoted dil. here and assumed to exclude boundaries at infinity), so that V,' is an open set. 
The concept of a dual access kernel is then easy to define: it is a kernel with two arguments where one 
argument is in Q and the other is in fl': Kd = Koix, q) where x £ 51 and q G f2'. A common example of a 
DAK arises in conjunction with the integral form of Poisson's equation for R", where the volume integral (or 
area integral in R^) of a density function times a DAK integrated over a source region produces the potential 
at some point outside the source region. When a DAK satisfies a replication property it is called dual-access 
collocation kernel (DACK) and the associated kernel space and inner product structure are known as a 
DACKS. Clearly a standard reproducing kernel, such as a Bergman kernel, is not an example of a DACK, 
since both the kernel arguments of a DACK are not even in the same subregion, but in general the occurrence 
of a DACK is tied to that of some reproducing kernel. A DACKS generalization of the reproducing property 
exists and is known as a replication or generalized least squares collocation property. For all currently known 
DACKs this replication condition can be expressed through the action of an associated one-to-one "involution 
mapping", p{q), that maps the region fl' into the interior of fi: p{q) G V.\ dft for q G ft'. It also satisfies 
the property p{p{q)) ~ q. The point p is called the involution point. For q g C this involution mapping 
becomes a "conjugate involution mapping" denoted p*{q) where p* is called the conjugate involution point. 
It may seem circular to specify Kd{x, q) with q G fl' and then apply a mapping of the form p{q) G ft so that 
functions are finally evaluated in 11; however, this approach has certain innate advantages since it allows, 
among other things, for the natural consideration of inverse source problems. As an example, consider the 
first DACK studied in the sequel; l/{z — Zk), where Zk is inside the complex unit disk and z in outside 
this disk and the associated conjugate involution is point is 1/z'l. Obviously this kernel form links a given 
(simple) pole region or source region with an exterior region of interest — i.e., a field region. In general, it is 



• kernel space 
• RKHS 



• DAK space 

• DACKS 



• Bergman kernel space 

• Szego kernel space 



• DIDACKS 

• SIDACKS 



10 



useful to limit the source region (i.e., the range space of here) to be a fixed proper closed subset of Vl' so 
that all the kernel forms studied are bounded. 

There are some minor complications here compared to RKHS theory. DACKS norms and kernels are no 
longer in general one-to-one correspondence with each other as they are in RKHS theory and a particular 
DACKS norm may have many different point replication kernels associated with it. Likewise for a single 
fixed kernel, DACK replication properties may hold for several different norms that are all specified over 
the same given (field) region of space. For example, as noted earlier, it will be shown that l/(z — z^) has 
two separate spaces or norms associated with it and the field region for both of these norms is the same. 
This non-uniqueness causes some nomenclatural problems since it means that the standard RKHS naming 
convention mentioned previously cannot be followed. This issue is circumvented here by simply specifying a 
norm structure along with the particular DACK of interest (which is generally understood from the context). 
As one might gather from the discussion given in Section [1] many questions associated with overall DAK 
and DACK structures are largely open and one goal of this series of articles is to attempt to motivate others 
to explore these structures. Much of the nomenclature just outlined is aimed at incrementally advancing 
this implicit goal under the assumption that labeling is a useful first step. 

Given the potential connections between RKHSs and DACKSs outlined above, it is natural to wonder 
about the specific connections between well known reproducing kernels and the particular DACKs examples 
studied here. At the end of Section [5] a general connection will be shown to exist between the primary 
replication property for the first norm studied and Cauchy's integral formula. It is well known that the 
reproducing property of Bergman and Szego kernels for the interior of a disk can be also be tied to Cauchy's 
integral formula. Along similar lines it is fairly easy to show that for the exterior of a unit disk the primary 
DIDACKS replication property is linked to Cauchy's integral formula, as is the Bergman kernel reproducing 
property for the for the same region. Obviously an indirect link exists between all these kernels. Other 
direct connections are established in the sequel between DACKS kernels and the reproducing kernels just 
mentioned. For the exterior problem at the end of Section [5] a direct link is shown to the interior Szego 
kernel. At the end of Section [6] it is shown that for the exterior of a unit disk the DIDACKS formalism can 
be linked to exterior Bergman space structures. Finally, from the material in Section [9] a direct connection 
between the interior form of the first norm and Szego kernel theory can be shown to exist. 

3 Notational and Allied Issues 

This section addresses various notational issues and summarizes various acronyms and terms used elsewhere 
in the paper. It is broken into various natural groupings where each item of interest is broken out separately. 
Often the motivation behind a particular notation choice is given since choice of notation is important, but 
rarely discussed. 

The two conventions discussed first are not standard, but are adopted here since they have strong edu- 
cational overtones. These educational overtones that are briefly considered in context: 

• For notational and conceptual clarity it is useful to be able to distinguish between the two normal uses 
of complex conjugation: (1) As an operator. (2) To denote independent complex conjugate variables. 
This distinction can be easily made since there are also two symbols in common usage: the superscript 
* and overbar. Because the matrix transpose operator is usually denoted by a superscript T, which 
is the convention adopted here, by analogy it is natural to use a superscript * to denote the complex 
conjugate operator (i.e., /*). From a matrix perspective this usage is consistent and allows a dagger 
superscript to denote the Hermitian conjugate of a matrix; i.e., for Hermitian matrix H it follows that 
ijt {^H^y ^ ^ This superscript * is the default usage here, while an overbar is reserved 

for those rare occasions where there is a real need to introduce independent conjugate variables (i.e. 
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z and z). Most remaining ambiguities can be removed by the judicious insertion of parentheses [for 
example, see the right hand side of (dH)] where the issue of the independence of z and z* does not arise 
(nor even need it be considered). Further, any residual confusion that remains in specific functional 
forms can be eliminated simply by referring to the series form that is used to define the admissible 
class of functions ([9]). 

Conjugate variables are useful and can occasionally greatly simplify an argument. They are also widely 
used in the analytic function literature (such as in conjunction with Hardy and Bergman spaces). The 
issue of the legitimacy of conjugate variables was brought up and settled long ago [20], but the concept 
of conjugate variables is not generally covered in introductory books on complex variables. By using 
both symbols in the way each is most naturally used, the full advantages of both notations can be 
accrued. 

Another approach often taken in the literature is to introduce the Wirtinger differential operators [26] 
commonly written as 



Other labels associated with these operators are Wirtinger derivatives and the Wirtinger calculus. 
Here the concept of an antianalytic function can be introduced to restore the symmetry between 
complex and complex conjugate variables and thus to fully leverage the innate advantage of these 
Wirtinger operators |12j . but this approach is not taken here. (The "mirror symmetry" of analytic 
and antianalytic function theory has a bearing on the next item under discussion, but this possibility 
is also not followed up in the sequel.) 

The reader who is an expert in analytic function theory may not see the utility of two separate 
conjugation symbols since the underpinning conceptual issues and tools are so obvious to him or her, 
but such a reader might want to consider the following hypothetical examples from matrix theory that 
suggests the sort of issues that a neophyte may encounter upon first exposure to analytic function 
theory when an overbar is used as an operator. Specifically, suppose that an overbar is used in matrix 
theory to indicate transposition and that matrices that are functions of other matrices are considered: 
M(N), for matrices M and N. Then what is the exact meaning, if any, of M(N), of M(N), of M(N), 

of M(N), of M(N), of M(N), of M(N)-\ of M^^N) or finally of NFT(N). Moreover, what is the 
exact difference, if any, of each of these expressions and can the reader easily interpret a handwritten 
version of these expressions? That any effort whatsoever is required to answer these sort of questions 
indicates a possible inherent ambiguity in the notation itself. 

Finally, the use of two separate conjugation symbols means that functions or expressions that are not 
analytic can be manipulated and handled just as naturally as ones that are. Specifically, the resulting 
expressions for inner product evaluations need not necessarily be analytic. 

• For pedagogical reasons, it is useful to distinguish between the two conventions that can be used to 
apply complex conjugation to the arguments of a complex valued inner product. An inner product (or 
norm) where conjugation is applied to the first argument will said to be in left-conjugate form (LCF), 
while if conjugation is applied to the second factor it will said to be in right-conjugate form (RCF). 
As an example, consider the complex form of the real inner product that was introduced in Section [2] 
to explain reproducing kernels, where T := [a, b] is now taken to denote a line segment along the real 
axis of the complex plane. The LCF of this inner product is (/, g)LCF '■= Jxifi^))* '^^j while the 
RCF is (/, g)RCF ■— Jjfiz) {g{z))* dx. Generally, when deahng with "analytic functional analysis" 
(i.e., complex RKHS theory, Bergman kernel theory, Szego kernel theory, theory etc.) it is the 




and 




(1) 
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custom to use RCF inner product structures, which is the convention originally adopted by Bergman. 
The LCF is adopted here for several specific reasons: 

— Here the resulting expressions from an LCF inner product application may not be analytic, but 
this does not matter in an applied setting. 

— The linear equation sets that result from an LCF treatment are much more intuitive and easily 
derived [c.f., jS])]. 

— LCF is the standard convention in quantum mechanics and elementary linear algebra. 

— LCF has also occasionally been adopted by mathematicians tlSj or used in concert with RCF by 
mathematicians 0]. 

— From an educational perspective, the LCF convention allows much of the intuition gained in matrix 
theory to be directly leveraged. Later, after the student is more comfortable with the concepts 
involved, the RCF can be introduced. (It is also worth noting that in the past, compositions 
in the reverse order have been periodically tried in introductory group theory courses and other 
areas and have been generally found to be more confusing to students than the standard operator 
ordering, which they encounter during a first exposure to calculus and matrix theory.) When a 
reproducing kernel, K^^ is considered as an operator it is natural to place it on the left in analogy 
with matrix theory; moreover, in the complex setting, at a certain level, it is largely a matter 
of choice as to whether the end result of [Kji, /) or (/, Kb) is required to be analytic (c.f, the 
concept of antianalytic functions mentioned above) . 

— From an applied perspective the choice of LCF or RCF is strictly a matter of convenience since 
either choice can be justified. Since matters are treated somewhat explicitly here, no inconvenience 
or confusion results by using LCF. In either case some care is always required in using the resulting 
expressions (for example, see the above discussion of independent conjugate variables). 

— Obviously, it is easy to convert from LCF to RCF since (g, f )RCF {g, IYlcf- Further, while 
the RCF may be a real conceptual imposition to the neophyte or to the practitioner from another 
field, the conversion from one form to the other should pose no real problem for professional 
mathematicians — especially when one form or the other is maintained throughout a particular 
sub-discipline (as mentioned above, it is a working assumption here that DIDACKS theory can 
be treated as an independent sub-discipline). 

With this notation, the real form of an inner product is given by ^(g, f)RCF + ^{g, f)LCF- The 
paper's formalism was tested and checked in LCF. In particular, looking ahead notice that retracing 
the derivation leading up to ^ employing RCF yields fi^Tj^f^p = aJ^^^, where Trcf denotes 
the matrix whose elements are Trcf^ ^ '■= {Bk' , Bk)RCF and Arcf the vector whose element are 
^RCFk [Bk, Drcf — which is equivalent to, but slightly more cumbersome than, ([5]) [and thus 
possibly confusing to a novice]. 



In the sequel the labels real and complex setting are used in the following context: 

• The real setting the setting where the function to be approximated and the approximating 
function, as well as the basis functions that compose it, are real valued and have arguments in M". 
The inner product is also real valued so that any scalar quantities involved are also real valued. 

• The complex setting the setting where the function to be approximated and the approximating 
function, as well as the basis functions that compose it, are complex valued and have arguments in C. 
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The inner product is also complex valued and unless otherwise stated it is assumed to be in LCF. It is 
also assumed that scalar quantities are also generally complex valued. 



Before summarizing various acronyms, it is useful to discuss some of the index and subscript conventions: 

• The subscripts k and k' range over the values 1, 2, 3, . . . , A^fe . 
The subscript conventions associated with the major norms studied are: 

• A subscript a is used to denote the standard norm where a =^ the unit circle = the boundary of 
unit the disk (9f2) = the inner boundary of the exterior of the unit disk. (While the outward pointing 
normal for the interior and the exterior of the unit disk differ in sign, this normal is not explicitly 
introduced here for any of the standard- norm usages so this equivalence strictly holds.) 

• Since the inner products associated with the Dirichlet integral assume different forms in the complex 
and real setting two different subscripts are used to distinguish these settings: 

— D =^ complex setting 

- E real setting. 



There are also a number of standard acronyms worth summarizing here: 

• LLSQ = linear least squares (when no confusion can arise this acronym is also used to denote a 
generalized linear least squares setting) 

• NLLSQ = non-linear least squares 

• RKHS = reproducing kernel Hilbert space 

• SVD = singular value decomposition. 

Finally it is useful to summarize the main acronyms that are introduced in the paper: 

• DACK = dual-access collocation kernel 

• DACKS = dual-access collocation-kernel space 

• DAK = dual-access kernel 

• DIDACKS = Dirichlet integral DACKS (in either the real or complex setting) 

• LCF = left- conjugate form 

• RCF = right-conjugate form 

• SIDACKS = standard integral DACKS (in either the real or complex setting). 

Finally, as a linguistic aside, it is assumed herein that an acronym denotes an exact word sequence and 
thus, for example, that where it is appropriate the plural of an acronym can be formed by simply affixing 
a small letter s to the end of the acronym. Also, although there is no established uniform convention in 
the literature with regards to forming the plural mathematical symbols, here and in future articles in the 
sequence an apostrophe s is invariably used in order to provide a useful visual clue as to where the math 
symbol itself ends (i.e, a^'s is prefered over a^s). 
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4 Foundations of DIDACKS Theory 



In what follows the basic strategy is the usual one of approximation theory, which involves introducing a 
norm and then minimizing the cost function 

^■.= \\f-^r = \\fr-2{f,^) + M\ (2) 

where if is the approximating function and / is the function to be approximated. When ip is assumed to be 
a linear combination of fixed basis functions of the form 



(p{z) ^^HkBk{z) (3) 



k=l 

a (generalized) LLSQ problem results, which is the main focus of this article. 

For the complex setting z, iik and Bk G C and the inner product is assumed to be in LCF. For the real 
setting z e R" and f{z), Bk{z) and /x^ real valued. 

For the real setting, the values of ^Xk that minimize $ = ^{^jLk) can be found simply by solving the linear 
equation set that results from setting the partials of $ with respect to ^k' to zero (and diving by two): 

J2iBk',Bk)tik^iBk', f) , (4) 
fe=i 

which can be written more compactly as 



T/i = A (5) 



where T denotes the matrix whose elements are Tk'^k {Bk', Bk), A the vector whose element are Ak := 
{Bk, f) and /i the vector whose elements are fik- Equation ([5]) will prove to be pivotal in what follows and 
it occurs in different contexts at several places in the sequel. [Here for the complex setting it is a simple 
matter to show that minimization of $ = ^{f^k) again leads to the same expression — see the steps leading 
up to ^ or (|T7l).] 

The choice of inner product structure and of basis functions are of fundamental importance. As mentioned 
in Section[l]the first class of basis function studied here will be Bk{z) = l/{z — Zk), which is a DACKS kernel 
in the complex setting. In the real setting this basis function corresponds to a dipole term. In Section [T] 
it was also as noted that logarithmic basis functions are DACKS kernels in both the complex and real 
setting. 

Next, consider the Dirichlet integral and its implications here. Although the approach in the core part 
of the article is to consider the complex setting first and then transition to the real setting, in order to gain 
a complete perspective the reverse course is taken in this section, hence the setting will be considered 
exclusively before transitioning to the complex setting. Here a distinction is made between the Dirichlet 
integral and Dirichlet form — although it is not the usual custom in the literature. For the real setting, the 
term "Dirichlet integral", denoted D[/, g], is taken to be an integral over a class of admissible functions that 
has a mathematical form proportional to the field energy (i.e., energy density) over the region of interest; i.e., 
([55)1 for the exterior of a unit disk in . For exterior disk regions the class of admissible functions here are 
essentially assumed to be bounded harmonic functions (with bounded partials of all orders :— C^) that tail 
off sufficiently fast at infinity. When a function F , which is not assumed to be harmonic, satisfies much less 
restrictive conditions, the well known Dirichlet principle can be considered and the term "Dirichlet form" 
will be reserved for the associated form D[i^, F] in this context. Looking ahead to the complex setting, it 
is important to keep this distinction in mind since the complex form of Dirichlet 's integral adopted here is 
derived under the explicit assumption that all admissible functions are analytic and, in this sense, it is not 
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a Dirichlet form, nor is the question of the proper definition of a Dirichlet form in the complex setting even 
raised in the main body of the text. (Ahhough this is an intrinsically interesting question and it is worth 
observing that at the genesis of complex variable theory properties of analytic functions and other results 
were obtained using Dirichlet forms through the action of Dirichlet's principle.) 

Next, recall that the Dirichlet principle basically asserts that the minimization of the form D[F, F] over 
some specified interior region, subject to the constraint that F match given Dirichlet boundary conditions, 
implies F is harmonic. Here for Dirichlet's principle to hold F must satisfy certain minimal conditions (i.e., 
least restrictive conditions) , as must the specified boundary conditions and the shape of the boundary itself. 
Dirichlet's principle is also assumed to hold for exterior regions, but as discussed below there are a few more 
complications so this case is usually not addressed in the broader literature. The historical counterexamples 
[5] for interior regions, their associated history [TB], the general nature of the requirements [131 [22] and typical 
recent associated research efforts |13) are of great interest in their own right, but are they are not especially 
relevant here since Dirichlet's principle is a tautology in DIDACKS theory — harmonicity is already assumed 
in the class of admissible functions itself. Clearly there is a middle ground between these two stances where 
harmonicity is not assumed for the class of admissible functions, but the assumptions are strong enough 
on the class of admissible functions and the regions of interest so that Dirichlet's principle can be derived 
simply without untoward complications. Historically assuming a certain class of explicitly defined admissible 
functions has been considered the acceptable way to attack Dirichlet's principle, with the first use in this 
context going back to Dirichlet himself [T51 p. 40]. Since the focus has generally been on proving Dirichlet's 
principle under the most general assumptions, these original attempts were flawed and counterexamples 
played a major role in uncovering the conceptual errors that were made at various historical junctures. If the 
boundary of the domain under consideration is not pathological, many of the innate difficulties in a proof 
of Dirichlet's principle arise from the broad class of allowed functions that are inherent in the Hilbert space 
concept — especially completeness with regards to infinite sequences of functions. Assuming a well formulated 
structured pre-Hilbert space setting by-passes these problems as well as problems in transforming Dirichlet's 
integral into an acceptable inner product structure. For a structured pre-Hilbert space setting, it is first 
of all assumed that a pre-Hilbert space setting holds — i.e., one where only an inner product structure is 
assumed — and that a "structure" is imposed on it by the requirement that all functions be drawn from some 
specified class of admissible functions. Finally it is worth noting that, as defined here the energy norm and 
Dirichlet integral are equivalent to within a factor of proportionality. The value of this constant factor is 
unimportant for minimization problems so long as it is introduced and handled consistently. 

Given that currently a pre-Hilbert space setting is not generally the first-line choice, its adoption here 
is worth discussing and justifying at some length. There are three primary reasons that a structured pre- 
Hilbert space setting is used here: (1) ease of use issues, (2) non-uniqueness issues and (3) mathematical 
consistency issues for energy based inner product structures that are defined over unbounded domains. The 
ease of use issues can be readily dispatched with: The main point is that various mathematical derivations 
can be readily carried out from the given properties of the class of admissible functions. Along similar lines 
it should be evident that the pre-Hilbert space setting employed here has certain innate educational benefits. 
(In this regard [T^ gives a fairly complete introduction to RKHS theory based on a pre-Hilbert space setting, 
but at a higher level than assumed here.) Each of the remaining two issues will now be addressed in turn. 

With regards to the second point, in the real setting is well known that various continuous source 
distributions inside some bounded region O' produce identical fields in the complementary unbounded region 
VI and thus one might expect ([5]) to be not only ill-conditioned, but generally ill-posed. What is less well 
known is that a proof exists in the geophysical literature that point source basis functions are independent — 
or, in other words, that a localized finite number of non-zero bounded scalar point sources inside some given 
fixed sphere of M", for n > 2, cannot exist that produces a null field in the exterior of the same sphere |24) . 
( [25] discusses a different type of inverse source problem for logarithmic potentials where the shape of 
the source region is to be determined.) Obviously this implies that ([5]) is theoretically invertible for point 
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source basis functions. 

A full discussion of point source uniqueness will be given in another venue, but observe that by using 
through (fT51) it is a simple matter to deduce that point pole distributions in both C and are unique. 
In particular, consider the kernel l/{z — z^) in the complex setting first and suppose, to the contrary, that a 
pole distribution of the form ip — X^fc^i ^i-kjiz — Zk) exists where Lp{z) = holds for all z e 51 (i.e., \z\ > 1) 
and all the source terms are assumed to be nonzero. Thus ^tT/i — must clearly hold for any 

choice of norm. Since T is a (Hermitian) Gram matrix all of its eigenvalues are real and nonnegative. Let 
Afc denote each of these Nk eigenvalues, then |T| = n£!i > holds for any norm. Combining this 
result with the fact that p9)) implies |T| 7^ holds for the standard integral norm setting, it is obvious that 
Afc > and thus that T is positive definite in this same norm setting. Then, by either expanding /i in terms 
of unit eigenvectors and substituting the result into fj}T fi, or by diagonalizing T in the same form, it is 
apparent that || </? ||cr = /z = . This contradiction proves the desired general assertion in the complex 
setting. For the real setting the uniqueness proof for dipoles is even easier since it follows immediately 
from the above result when coupled with the facts that Re {l/(z — z^)} — {x — Xk)/ {x — Xk)^ + (?/ ^ Vk)"^, 
Im {1/ {z — Zk)} = {y — Vk)l \/ (2; — + (2/ ^ VkY ^ii'i that both the real and imaginary parts of a complex 
function are harmonic. 

While it may be true that point source or point pole distributions are unique, ill-posedness considerations 
clearly remain and as a practical matter these theoretical existence and uniqueness proofs do not in any 
sense guarantee viable solutions. Thus theoretical non-uniqueness issues are by-passed here and a utilitarian 
perspective is adopted. In this perspective, the minimal requirement is that be (accurately) numerically 
computable, which is the assumption made here. This immediately shifts the perspective from abstract 
considerations to condition number and other implementation concerns. As such, condition number issues 
will be addressed more fully later at the end of this section. Also this and other implementation concerns 
are addressed in Section [51 To see the inherent difficulties here in a Hilbert space setting consider the limit 
of a sequence of approximating functions given by \im.N^~frxi Vn^. i where is specified by the right hand 
side of ([3]) with a set of fundamental solution basis functions for Bk- Clearly, in this limit it is possible 
to approach continuous distributions that are known to be indeterminate. This is one of many troubling 
sequences whose limits must be admitted from the completeness requirements placed on a Hilbert space and 
which a pre-Hilbert space setting by-passes. 

The third point is that various conditions must be enforced to transform a Dirichlet integral into a proper 
inner-product structure and a structured pre-Hilbert space setting is the natural way to do this. Historically 
speaking, the assumption of a structured pre-Hilbert space setting has been very common in dealing with 
energy norms, especially for unbounded regions |;5j. There are four primary problems to be overcome. First, 
the functions to be approximated and the approximating functions must be well behaved on the boundary of 
fi. Second, D[i^, F] is inherently indeterminate; for example, in the real setting consider the fact that W{F + 
C) = VF for any constant C [where V denotes the usual gradients :— {d /dx, d /dy)'^]. This problem 
is easily solved by imposing the requirement on the class of admissible functions that limi^i^oo Fi^) = 0. 
Assuming that F and G are sufficiently smooth, the third problem is that D[i^, G] many not be bounded. 
Assuming that WF and VG are bounded, this is easily solved by requiring that WF and VG also tail-off 
fast enough as distances approach infinity. Fourth, various smoothness issues are present that can easily 
be by-passed by imposing harmonicity or by requiring analyticity in the complex setting. As discussed 
below, since it is fairly easy to solve all of these problems in the complex setting by employing a power 
series representation, this is the approach taken here. (Researchers can simply infer the appropriate minimal 
requirements, while the student reader who wishes to consider the required properties from a more concise 
perspective should consult a standard reference such as [2] and lecturers should see Item 9 of Appendix A, 
where the motivation for the approach taken here with regards to such matters is presented.) 

Three other observations are relevant. First, it is also easy to characterize the class of admissible functions 
by assuming that they are all generated by well behaved density functions. Second, for bounded domains it 
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is common to employ a Sobolev space instead of a structured pre-Hilbert space to overcome these Dirichlet 
integral issues, but this entails a basic modification of the inner product structure itself; i.e., in the first 
Sobolev norm setting a volume integral over ^l of the product of F and G is added to D[F, G] and it is then 
required that the resulting inner product evaluations be bounded. Third, various properties defining the 
class of admissible functions need not be concisely stated and, in fact, it may be convenient to have a certain 
amount of redundancy built into them. (As elsewhere in this article the reader is reminded that generality 
and conciseness are not to be confused with rigor, as some are prone to do, since rigor does not necessarily 
require either of these other two conditions.) The point is that if ease of use (or exposition) is adopted as 
a primary goal, then one property or another can be more effectively employed to perform some particular 
derivation so long as it is clear that the functions involved are inside the overall common intersection of 
subspaces that results when all of the properties concurrently hold (generally each property allows linear 
superpositions so that a finite span of functions results since linear combinations of functions that satisfy a 
property, also satisfy the property). 

Clearly certain types of functions to be fit may fail to fall within the parameters of the given class of 
allowed admissible functions. These exceptional cases can generally be handled by applying special attention 
to them and introducing modifications on a case by case basis as needed. It is instructive to examine two 
concrete examples along these lines. Towards that end, consider the complex setting, which is the main 
one assumed in the sequel. Here the properties for the admissible class of functions are specified by ([9|). 
but it is clear that there are common analytic functions over the region of interest that do not fall within 
the requirements of this class of admissible functions. First, although /(z) = In z is analytic in fi, it fails 
to meet either of the criteria specified in ([9]). In Section [7] it will be shown that both these difficulties can 
be overcome by introducing the idea of a paired point. This paired point formalism physically amounts to 
considering "finite extended poles" only. As a second example, consider f{z) = sin z, which will be analyzed 
in more detail in Section [5] and it will be shown there that this function can be modified to produce a 
suitable function. This general issue of admissible functions will be dealt with in more detail in Section [5] 
Here it is worth emphasizing again that in the sequel the class of admissible real functions is obtained from 
the class of analytic functions ^ by simply taking the real part of this series and as discussed elsewhere 
this correspondence is one-to-one since such admissible real functions (aside from branch-cut related issues 
for certain cases) have ([9|) as a unique harmonic completion. This correspondence is essentially an easy 
way of circumventing the various issues mentioned above in conjunction with Dirichlet integral based inner 
products. 

As indicated in the introduction and in Appendix A, it is assumed that all readers have at least some 
understanding of functional analysis, but for those readers who are not otherwise comfortable with Dirichlet 
forms or integrals in a few additional comments are perhaps in order. First the actual series representation 
of the class of admissible real valued functions for the exterior of a unit disk, which are denoted here by 
F{x, y), can easily be explicitly determined as follows. In what follows x and y are Cartesian coordinates 
and r and are the corresponding polar coordinates. Each term of ^ can be rewritten as 

a„ a„ e^*"® a„ cos n9 — iun sin nO 

for n > and \z\ > 1, with z G C, a„ €E C and n an integer. It thus follows immediately by taking the real 
part of ^ and reidentifying coefficients that 

F:.Re{/(z)}:= Re { £ a^z^ } ^ g + ^" "V (6) 

n— 1 n— 1 

where a^, 5^ S M. Observe that © is just a standard harmonic series representation over the exterior of a 
unit disk of a harmonic function that vanishes at infinity. One important property of any admissible F is 
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that it must be bounded with respect to the inner product and norm structure: WFW^ := (F, F) < oo. This 
is accomplished here through the properties of the series forms ([B]) and hence, some natural but not 
otherwise obvious constraints must also be placed on the coefficients in these series representations so that 
any coefficient series that is the result of an inner product or norm evaluation converges. 

Next consider the issue of explicitly showing that Dirichlet integrals of ([6]) are bounded over the 
exterior region of interest: := { x, y \ r — \/ + ^ 1 }■ I3efore considering inner products based on 
Dirichlet integrals consider the following inner product, which will be introduced at the end of Section [5] in 
a different setting: 

{F,G)^^ II FGdA, (7) 
n 

where dA is the element of area, so that dA — dx dy in Cartesian coordinates and d A — r dr dO in standard 
polar coordinates. This inner product, however, is unsatisfactory for the class of functions considered here 
since even the first term in the series on the right hand side of ([6]) leads to an unbounded integral: 

oo 27r 2 

||sin0/r||5j= I I ^^rdrde (8) 

r=l 9=0 

(which is clearly unbounded). Alternatively, since in polar coordinates 

dF ^ 1 dF ^ 
\/F = — + -TTTT ee and thus |Vi^p = 
or r oO 

it is clear that ([35]) will result in a bounded integral for well behaved a'^ and b'^ so that a consistent norm 
and inner product structure can be based on ^ for Dirichlet integrals. 

To guarantee boundedness of inner product and norm evaluations, it is also necessary to place restrictions 
on the series coefficients of ^ and ©. The actual properties that the coefficients a„ in © or a'^ and b'^ in 
dH) must satisfy can be easily found by direct substitution. For example, if (O is substituted directly into 
D[F, F] given by (l38|) then the resulting integrals can be directly evaluated by the same procedures used in 
Section m This is not done here to emphasize the point that so long as the associated norms are bounded 
not only are these actual expressions not required here, but that all that is required is that a series of the 
form (O or ([H]) must exist in principle so that even the form of the coefficients a„ in ^ or of a'^ and b'^ in 
(|6l) need not be known explicitly. 

Finally, from an applied perspective the requirement that solutions to ([5]) be numerically computable 
is clearly significant for both inverse source determination and functional approximation problems, but 
condition number considerations differ somewhat for these two areas of interest. In modeling problems 
one wishes to minimize the modeling error directly and by introducing a cost function that minimizes this 
modeling error it can be done effectively; moreover, the minimization of ([2]) that is implied by ([5]) ensures 
that the modeling match is good regardless of the size of the condition number that is associated with the 
underlying matrix equation set (at least up to a point). The underlying issues here can perhaps be most 
clearly understood by way of an analogy with a standard LLSQ process (i.e., not a generalized LLSQ one), 
where a Euclidean sum of squares of errors is minimized and the linear system of equations that results 
is, say Ax = b. Here the vector b is known and the ill-conditioned matrix A is known, but the vector of 
parameters x is unknown and is to be determined. Let x denote the true value of x and x the numerical 
solution. Then as an optimization criteria, clearly there is a large difference between requiring that |Ax — bp 
be small for a successful solution and requiring that |x — be small: it is this difference that underlies the 
point made above. The underpinning issues here will be addressed further in Section [8l as well as in other 
articles in the sequence. 



dF 


z 

^ 1 


dF 


dr 
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5 Standard Integral Norm Setting and General Motivation 



This section delineates the basic SIDACKS approach. In what follows it is assumed that functions are 
bounded and analytic in the exterior and on the boundary of some unit disk specified by f2 := {z G C | \z\> 
1}. It is further assumed that the magnitude of such functions fall off at least as fast as l/\z\ as \z\ oo. 
(As mentioned in Section [4] this condition is somewhat restrictive — for example, if lim|2|_^oo f{z) — constant, 
then the constant term can be subtracted off and later added back on after the fit has been performed; 
however, even after this modification many analytic functions will still fail to meet this general criteria [see 
Section [5] for various examples].) For \z\ > 1, such functions can be efficiently represented by a power series 
of the form: 

oo 

/(z) = ^a„z-" =^ lim /(z)^0, (9) 

where a„ e C. It is assumed that the series in ^ is well behaved (i.e., that it yields finite norms) and that 
it can be operated on termwise to yield series that are bounded and point- wise convergent, which is to say 
that 11/11 < oo and that the power series for F{z) = f{l/z) has a radius of convergence that is less than one. 
Equation ^ essentially specifies the admissible class of analytic functions that will be considered. It has 
been stated in a somewhat redundant fashion for emphasis (as discussed in Section|3]and Item 9 of Appendix 
A, redundancy here is not necessarily considered undesirable). In the sequel, although this power series is 
assumed to exist, none of the actual terms in the series are explicitly required provided that there are no 
poles located at the origin. While ^ has been written in terms of the function to be approximated, the 
approximating function is composed of linear combinations of basis functions, which are simple poles with 
specified positions that clearly satisfy jl]): 



^ = E^- (10) 



Here it is assumed that iik ^ 'C and that < \zk\ < 1. While Zk 7^ is assumed for convenience, this 
restrictions can easily be removed (see Section [8]). Using the geometric series, each of the pole terms 



appearing in (jlOp can be reexpressed as 



^ (11) 



Zfe ^ z"+i 
n=0 



As discussed in Sectional the basic strategy is to introduce an appropriate LCF norm structure and then 
to minimize $ := ||/ — — — 2(/, ip) + ||v5|p. Towards that end consider the following inner product: 



2ir 

(/,5). / \r^^)9^^)\.J^ ^ (12) 

9=0 

where polar coordinates have been introduced (i.e., z — re*^, Zfc = rfcC*^*). [The factor of l/(27r) has been 
inserted here since it is clearly inappropriate to lump it in with the kernel l/(z — z^) — as is the usual 
convention in RKHS theory — and the alternative is to allow factors of tt in the point replication property 
itself.] Clearly when |1/ - := (/ _ / _ <^)2 = (|/ _ _ (p|)2 || |j _ ^| ||2 minimized 

the boundary conditions on the unit disk are also minimized and thus f ~ f follows immediately from the 
uniqueness of Dirichlet boundary conditions. Replacing 5 by l/(z — z^) in (fT2|) and using ([9]) and (fTTjl yields: 

27r 

„ 00 00 00 

{f,{^-^k)-')„ = ^ / E z]:'e'^^-'+'^ue^-Y,<4 ■ (13) 
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Introducing pk — 1/zk allows (jl3p to be rewritten in the final form 

{fAz-zkr'),=pk{f{pi)r ■ (14) 

When (fTO|) is used to reexpress cp and the result is substituted into <i>cr := ||/ — iy9||^, a bi-linear form 
in terms of Hk* and fik results, which can be further simplified by introducing {{z — Zk)~^, f)a- ■= Ak and 
{{z— Zki)~^ , {z— Zk)~^)c! '■— Tk'.k- The resulting expression for $0- can be further simplified since Ak and Tk\ k 
can be rewritten as closed form expressions by using (|14p . [Here Tk'.k = {PkP*k') I {PkPk' ^ 1) ^ (1 ~ ^fc-^fc')^^ 
and Ak = P*kf{P*k)-\ Thus it is only necessary to take partials of this bi-linear form with respect to the 
source parameters and then set the result to zero in order to obtain a closed-form equation set which gives 
the values of the pole strengths (/ij ) that minimize <i>cr- Specifically, if T denotes the matrix whose elements 
are Tk'. k j M the vector whose elements are /i^ and A the vector whose element are Ak then, as shown in 
Section [HI using complex conjugate variables, the linear equation set governing the source strengths can be 
written as 

T[„] = A[„]. (15) 

In p5)) a generally omitted [optional] a subscript is indicated: i.e., T[„] :— T, := ^ and A[„] := A. 

Equation (fT5)) can also be obtained without resorting to complex conjugate coordinates by following a 
more involved sequence of steps. First introduce variables for the real and imaginary components of Hk' 
ak '■= Re{/ifc} and fik '■— Im{/ife}. This yields a cost function $0- = '^a{otj, Pj) with real arguments. Next 
take partials with respect to aj and jSj and set the results to zero, which yields two real coupled linear 
equation sets that can be rewritten as a single real equation set in block matrix form: 

(16) 



Tr 


Ti 


-Ti 


Tr 




In ([16)) the subscript R on a vector or matrix denotes the real part of the vector or matrix in question, and 
I denotes the imaginary part. Also a denotes the Nk vector whose components are ak and f3 the Nk vector 
whose components are (3k- Finally, by substituting T = Tr +iTi, fi = a + i (3 and A = Aa + i Ai into 
it is a trivial matter to show it is equivalent to (fTBl) ^2T]. 

Here T can be rewritten as the product of a matrix L whose elements are Lk' , fc = 1/ {pk — ) and a 
diagonal matrix. Dp whose elements are j pj : 

T = LDp . (17) 

The determinant of L, |L|, is the well-known Cauchy determinant. Introducing the standard symbols a; :— pi 
and bj := —z* used in closed-form expressions of Cauchy's determinant, the associated matrix has elements 
of the form l/{ai + bj) and the closed- form expressions itself can be written as [H p. 268] 

l\iat-aj){bi-bj) n (Pt - Pj ) (^j - z* ) 

n ia,+bj) n ip^~z*) 

While (HH]) is frequently stated for and bj G R, since the standard proof is by induction and involves normal 
algebraic manipulation |9j it clearly holds for and bj G C as well. Combining |Dp| — P1P2P3 • • 'PNk with 
this closed-form expression gives: 



Nk 



UiP^-Pj)i^ 

i>'j 



J ^ ^i) Nk 



Nk 

n {p^ - z*) '==1 
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Clearly |T| 7^ when < |zfe| < 1 since Zk 7^ Zk' holds for k' ^ k. Thus T^^ exists, but as discussed in 
Sections H] and [5] some care is normally required in solving (|15p or p6p . 

Higher order poles can be easily incorporated into the fitting function ip since 

rf" ( 1 \ ^ (-1)"'"^! ^ f^—] (20) 

dz^\z-Zk) (Z-Zfe)'"+1 dz]^\z-Zk) ^ ' 

and when partials with respect to z^ are taken on both sides of (fT4|) analogous closed-form expressions are 
produced for higher order poles. For example, a general fitting function of the form 



k—l m—1 ^ 



can be easily handled, where the Vk, m's are the different pole strengths of various orders at distinct locations. 
(Here Nm{k) is the number of different types of poles at location z^.) Clearly when this form is substituted 
in place of (jlOp into (f>o- := ||/ — a set of closed- form equations result for these various higher order pole 
strengths. 

So far the development has focused on the exterior of a unit disk; nevertheless, from an applications 
perspective a number of generalizations are possible. For example, if the region of interest consists of the 
exterior of a disk of different size that is centered over a different point then obviously a translation and 
rescaling can be applied beforehand and then, after a fit has been performed, the results can be rescaled and 
translated back to the original coordinates. Another fairly obvious approach is to employ (|15p to implement 
either a local or global interpolation scheme. Such schemes can be based on the fitting function forms given 
by pU]) or more general ones such as (HH). Obviously the specified data for such interpolation schemes — 
including the function itself and its derivatives up through order [Nm{k) — 1] — is matched at all of the 
conjugate involution points. 

As mentioned in Section [U the basic SIDACKS relationships — especially — are disarmingly simple 
and it is natural to wonder if and how they are related to known results. The most direct and obvious 
connection is of (|14p to the Cauchy integral formula. The Cauchy integral formula for the interior of a unit 
disk can be written as 

\w\ = l 

where |C| < 1 and F{() is analytic. Introducing z = 1/C, and relabeling F{Q := F{l/z) :— f{z) in ((22)) yields 
. Several RKHSs for either the interior or exterior region of the unit disk can also be tied to Cauchy's 
integral formula so the SIDACKS approach is either directly or indirectly related to them as well. Also 
as previously mentioned, the most significant RKHS connections are to the Bergman kernels for both the 
interior and exterior of a unit disk and to Szego kernel theory. Since Bergman kernels for are analogous 
to those for C, this later case is focused on in the sequel. These RKHS structures, of course, have RCF inner 
products and kernels. 

The link between the Bergman kernel for the interior of a unit disk (Kb) is worth briefiy considering 
since it too has links to the Cauchy integral formula. Let fl' := 17' U 90', then the Bergman inner product 
for the unit disk (0') is commonly written in RCF and is 

{F,G)^:^ II F{w){G{w)r dA (23) 



where the area integration is performed over the components of w with w G O'. The Bergman kernel is then 

Kb{C,w)^- \ (24) 

TT (1 — C,*WY 
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where, as above, C < 1. As will be shown next, this kernel has the reproducing property 



(F, Kb)^ = F(C) . (25) 
First observe that the derivative of the geometric series 1/(1 ~ x) — X!^o^^^ yields 



^ oo 



^ oo 

KB{(:,w) = -Y,{Cwr . (26) 



TT 

n=0 



C ^ w''K*B dA (27) 



so (1^^ can be rewritten as 



From (pS)) one can easily show 



by rewriting c? A in terms of polar coordinates and then performing the indicated integration over the unit 
disk, as was done, for example, in deriving (|13[) . Next, let 77 > 1, then substituting the geometric series for 
l/{r] — w) into ( • , Kb)^ and using ([ST]) , it follows immediately that 

(X^-u,)-\KB)^ = ^^y (28) 

To show the connection to the Cauchy integral formula consider a closed curve 7 which is completely outside 
the unit disk and suppose that F{z) is analytic inside and on 7, then the Cauchy integral formula for a point 
w which is inside the unit disk can be written as 

F{w) ^ ^ <f ^ drj . (29) 

1 

From ([28]) and (|29|) it follows immediately that 

{F,Kb)^= ^ iF{r,)({r,-w)-\KB) drj ^ ^ <f ^ dr^ ^ FiC) . (30) 
" 2tti J \ /B ZTTi J rj — Q 

7 7 

which is the desired reproducing property. Detailed rigorous discussions of Bergman kernel theory can be 
found in various entry level books [THl [5] • It is also worth noting that although Bergman spaces, which are 
an outgrowth of basic Bergman kernel theory, have been around for more than fifty years [4], they have 
continued to be an area of intermittent research that has fairly recently resurged [HI [11] . 

Just as was done with (j30|l for disk interiors, the Bergman kernel for the area over the exterior region of a 
unit disk can also be derived from the Cauchy integral formula. Thus the DACK for the standard norm and 
the Bergman kernel for the interior of a unit disk are linked through the Cauchy integral formula itself. It is 
also true that the Dirichlet integral based energy norm kernels, which are introduced in the sequel, can also 
be linked to the Cauchy integral formula, but in this case a more direct connection exists to the Bergman 
kernel for the exterior of a unit disk and this connection will be explored at the end of Section [Sj 

Links also exist between the standard DACKS formalism. Hardy space (H^) and Szego kernels. When 
equipped with the Szego kernel, can be considered a RKHS. In particular, the Szego kernel inserted 
into the standard Hardy inner product [denoted ( • , Oh] yields an alternative form of the Cauchy integral 
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formula; however, an even closer connection exists between the pre-Hilbert form of and the standard 
DACKS formalism. This connection is worth considering in more detail since the Cauchy integral formula 
need not be employed as an intermediary. 

First, in order to fix the notation it is useful to explicitly introduce the Szego kernel and its associated 
inner product spaces. Following Mate 16], let Hq denote the pre-Hilbert space of functions analytic in the 
interior of a unit disk (|2:| < 1) and continuous on the boundary \z\ = 1, with inner product specified by 



1 

27r7 



f{z)g*(z) (1/z) dz 



(31) 



|z| = l 



Here this line integral expression for the Hq inner product has been written in standard RCF to facilitate cross 
referencing. Since \z\ = 1 is understood in the integrand occuring in (|3ip . the relation dz 
can be used to yield the equivalent form 



I e 



ie 



dO = i z dO 



2t: 



27r J 

0=0 



(32) 



which is the standard form of the inner product for the pre-Hilbert space Lq over [0, 27r], with /(27r) — /(O) 
understood. The completion of Lq is, of course, the Hilbert space L^. Likewise the completion of Hq is 
denoted and it is an RKHS (given natural restrictions). As just noted the reproducing kernel for this 
RKHS is the Szego kernel Ks- Thus substituting g = Ks '■= 1/(1 — z s*) where s € C, with \s\ < 1, into 
(|32| produces 

(/, Ks)^2 - fis) . (33) 

The connections to standard SIDACKS theory are easily derived. While it is understood that SIDACKS 
theory holds for the exterior of the unit disk instead of the interior, this fact fails to show up explicitly in 
the inner products themselves (where |z| = 1 is understood), unless f or g depends on auxiliary parameters. 
Thus, as explained in the following, ((32)) is formally equivalent to the standard SIDACKS inner product given 
by mi), except for the use of RCF versus LCF. There are some subtleties here that involve the definition 
of admissible functions; however, these issues can be easily clarified after the fact. Equation (|55|) explicitly 
involves an auxiliary parameter (s) so some care is required. The procedure will be to recast (I14p from interior 
LCF into exterior RCF. The conversion to RCF simply entails taking the conjugate of the expression that 
results from combining (|12p and 



1 

2^ 



9=0 



(34) 



\z\ = l 



For clarity let z' :— 1/z, s — z^. and w := z' = r'e*^ = (l/^)e so that ([34]) can be rewritten as 



1 

2^ 



w*f{l/z')) 



'=0 



(1 — w*s) 



d0'^f{l/s)/s . 



(35) 



The function F{w) = f{l/w)/'w occuring on the right hand side of ([35]) is the relevant one and it is thus 
useful to consider the form of its power series in terms of w [which is based on ([5])]: 



n=0 



(36) 
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Since w* ~l/w for \w\ = 1, (j35p can be rewritten as 




L(l-u;*s)J i^i^i 



dff = F{s) 



(37) 



and this obviously corresponds to ([55]) . The subtlety indicated above in identifying the functions that occur 
in (fT2)) and ((32)) is thus clearly tied to the fact that / must actually be replaced by F. Notice that since 
l^l — 1 is understood in the integrand of (fT2|) . a factor of z*z can be inserted to aid in this reidentification 
process [the converse identification can easily be carried out using ()33p ]. 

Having shown a correspondence between inner product and replication kernel structures, the comparison 
between exterior SIDACKS theory and Szego kernel structures is now complete. 

For fits that do not involve derivatives (i.e., higher order poles) there is no real significant formal difference 
between SIDACKS fits using the simple-pole form l/{zk — z) and Szego kernel based interpolation. Thus for 
the interior of a unit disk it might be reasonable to consider the kernel l/{z}^ — z) as just the SIDACKS form 
of the Szego kernel; however, this nomenclature would be ambiguous here for various reasons that include 
labeling inconsistencies of norms and kernels — for example, at the end of Section [5] it will be shown that a 
different norm can be associated with a very similar kernel {\/[zk ~ z]" 1/^fc) over the same domain studied 
for \/{zk — z). The advantages of RKHS based approaches are well known, but even for interior problems 
DACKS theory is worth considering in its own right as an alternative to other interpolation schemes since 
it allows access to a number of diverse results in a easy, natural and seamless fashion. For example, for 
the SIDACKS interior case it is easy to use (j20|) to produce higher-order pole fits, but the same thing is 
cumbersome using the Szego or Bergman kernel. Also, as discussed in Section [U by providing a different 
perspective, consideration of aspects that might not otherwise be obvious can be entertained. Also as 
implied there, only a small portion of the overall DACKS formalism is contained in this article and much of 
the motivation for DACKS theory is contained in these other developments. Additionally, DACKS theory 
has a natural relationship to both energy minimization based approaches and to logarithmic source harmonic 
problems, which is generally missing from Bergman and Szego kernel theory since potential theory, inverse 
source theory and Dirichlet integral applications are not normally considered to be a part of Bergman or 
Szego kernel theory. 



As discussed in Section SI it is expedient to restrict G and H to be in some admissible class of functions so 
that D[G, H] can be associated with a positive definite inner product. Since the overall strategy employed 
here will be to recast ([55]) into a complex form that is more mathematically amenable, this class of admissible 
functions will be linked to this complex form. Thus the class of admissible functions given by dH) will be 
used to specify the functions of interest for this section and these functions are associated with the complex 
analog of ([55)1 and its various ramifications. After the developments of this section are complete, the real 
setting will be studied in the next section. The admissible class of functions introduced there will then 
be defined by the requirement that D[G, H] and its complex analog be linked through the act of standard 





(38) 



r=l e=o 
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completion of G and H (c.f., Definition [T]!.). The Dirichlet norm for the real setting is labeled the energy 
norm and it will not be actually introduced until late in the next section by (j66p . 

Thus consider the following known but underutilized complex generalization of D[G, H] over O C C: 



C30 2tz 



Dc[/, g] 



rdrdl 



(39) 



r=l 6=0 



where / and g are analytic in Q. Let df := f — ip = us{x, y) + i vs{x, y) define the modeling difference that 
is to be minimized, where us and vs are real valued functions of x and y. Since 5f is analytic, ug and vg 
satisfy the Cauchy-Riemann conditions and so 



dSfV /dSf 
dz ) \ dz 



dx 



dug 
dy 



dx 



dt's 
dy 



Thus minimizing 



Sf III) is equivalent to minimizing D[u, u] or D[u, v] , where 
||5/|l := V^cW,Sf] /27r. 



(40) 



(41) 



The associated inner product is, of course, given by (/, g) n :— Dc[f, g]/(27r). (Notice that $£i is real valued.) 

2 

The factor of l/(27r) here has been inserted for convenience and to help facilitate the identification of || 6f ||^ 
with field energy [which, aside from a units dependent constant, generally has a density of (Vi? • Vi/)/(27r)]. 
Taking the derivatives of ([9]) and (fTTjl and substituting the result into the right hand side of p9|) yields 



oo oo 



OO 27T 



(TO + l)a^^ie'("'+^^^ 



where m was reindexed to start at rather than 1. Clearly 

oo 

CXD „ ^ OO 

(/, (z-zfc)-i)^ = E(" + i)'<+i^^" / 7;:;rT5)i^^ = ^E"<^^" 

^■1—0 „„i n—1 



{n + l)z^e- 



-i(n+2)e 



„n+2 



■-( 

2zk V 



rdrde (42) 



df 



(43) 



Introducing ;= l/zk and fz{z) d f /d z(z), ()43p can be reexpressed as 



Pi 



fziPl) , 



(44) 



which is analogous to IH)) . As before, let Ak {{z — Zk) ^, /)d := (/, {z — Zk) ^Yq and T/ 
((z — Zfc/)"-'^, (z — Zfc)^-'^)^ , then $£) can be written as 



k'.k 



fc=ife'=i fc=i fc=i 

Equation (|44p can be used to obtain explicit expressions for and Tk'.k- 



-(Pk) fzipl) and Tk',k = 77 



1 (PfePfe')' 



2 (PfeP^ - 1)= 



(45) 



(46) 
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Introducing the matrix T and vectors A and /i, as in ^ or ([T5|) . along with the pair of complex conjugate 
variables fik and flk '■— fJ-k* into (j45p yields a cost function of the form = <I'd(m, m)- The minimum of 
is then found by setting d^o/dfij = (for j = 1, 2, 3 • • • , Nk): 

T,o, = A[o] , (47) 

in accord with (|15p . In (|47)) an [optional] 13 subscript is shown which is generally omitted: Tjdj :— T, 
:= /i and A[d] :— A. Obviously, applications here follow along the same lines as those of ()15|) so that, 
for example, fits based on ([20]) and (|21|) can be carried out here in the same fashion. 

For any suitable domain in the complex setting, when Dc [/, g] is defined analogously to ()39|) an intimate 
connection exists between Dc[f, g] and Bergman kernel theory for the same domain. For example, in analogy 
to (j23[) . when / and g are in the admissible class of functions, consider the following RCF Bergman inner 
product for the exterior of a unit disk: 

oo 27r 

(/,5)^:= I I J{z){g{z)rTdrde . (48) 

Since the derivatives of / and g always exist, introduce F{z) :— df/dz and G{z) := dg/dz. Then with 
Pk ■— l/zk and fz{z) :— df/dz this inner product yields 

(F, G);/2tt = (/, g)^ , (49) 

which clearly holds in general for any domain in the complex plane (with suitable restrictions on the class 
of admissible functions). 



7 Laplacian Inverse Source Theory in the Real Plane 

This section transcribes results obtained in Sections [5] and [6] for the complex setting to the real setting. 
A direct [23] , but slightly nonstandard, way to map a harmonic function specified over some region of 
M^, say W{x, y), into an analytic function fiv that is defined over the corresponding region of C is to set 
fw{z) Uw{x, y) + ivw{x, y), for real valued functions Uw and Vw , where 



Uw{x, y) 



dW 
dy 



and Vw{x, y) 



Using this correspondence in immediately yields 



dW 
dx 



(50) 



\dx ) \dy ) 



de 



(51) 



r=l 



so that minimizing — enforces Neumann boundary conditions on W. Note, however, that mini- 

mizing 1 1 /w I II) implies minimizing D[9W/9a;, dW/dx] or T)[dW/dy, dW/dy] rather than T)[W, W] itself, so 
that the analytic correspondence given by ([50]) will not be used in the sequel. 

As an alternative, thus consider the standard way to map W into an analytic function: 

Definition [TJl An analytic function fw{z) is said to be the standard completion of a harmonic func- 
tion VF, for some suitable common region in the x-y plane, when fw(z) := Uw{x, y) + ivw{x, y), where 
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Uw{x^ y) := W{x^ y) and is the analytic completion of Uw{x^ y) [i.e., Vw{x, y) is the unique complex 
function required to make fw analytic] . 

After taking branch cuts into account, the standard completion is unique since the real and complex parts 
of an analytic function are separately harmonic and there is thus a one-to-one correspondence of harmonic 
and analytic functions [5D] (for restrictions see Duren [TU]V 

For dipole and higher order poles in this standard completion process seamlessly meshes with the 
corresponding previous formalism given for pole strength determination in C and it allows for easy imple- 
mentations. For example, a dipole aligned along the a:-axis in is proportional to 

Q 1 X — Xk 

Uw '■— ———In — — — , (52) 

ox ^{x - XkY + {y- ykY ^/{x-XkY + {y-ykf 

which has a standard completion of 

, "'"^ + , '^'-''^ (53) 

V(x-a;fc)2 + (y-j/fc)2 ^(a;-a;fe)2 + (2/-yfe)2 z-Zk 

Likewise a dipole aligned along the y-axis in is proportional to 

d 1 

Uw ■= — T^ln — (54) 

dy ^(a:-Xkr + {y-ykr 

which has a standard completion of —i/{z — Zk). Since linear superposition holds for this standard completion 
process, an arbritrary combination of various dipole terms can be easily handled by the resulting DIDACKS 
formalism. Handling harmonic problems for the exterior of a disk thus presents no real conceptual 
difficulties since the usual procedure of mapping such problems into C and then recovering the results by an 
inverse transformation can be easily applied, or one can map the functions and inner products obtained in 
the complex setting into the real R^ setting by using relationships analogous to ([65|l and ((67|) . 

The inclusion of point logarithmic terms in R^ is not quite so trivial, but towards that end consider 
the R^ potential term Lk{x, y) In [l/-\/(x — Xk)"^ + (j/ — ^fe)^]- The standard completion of Lk{x, y) = 
ln[l/y(x^^a?^O^T~(y~-yfcF] = ln(l/|z-Zfe|) is 4(z) = In [1/(2; - z^)] = -ln(z-Zfc), which was not 
considered in Section [5] or Section [6l Nevertheless, retracing the steps that lead up to ((44|) for this logarithmic 
function is straightforward and yields the following suggestive relationship 

00 

(/, \n[{z-Zkr'])^ = := ^rbD ; (55) 



2 

n=l 



F 



where the symbol ~ indicates a formal manipulation (or relationship), which may or may not be fully 
justifiable. It is easy to see that ([55]) is not suitable as a framework for applications and then to discern 
the underlying reasons why. Towards that end, first substitute — ln(z — Zk') for f{z) in the left hand side 
of ([55|) (where k' ^ k). Next replace k by k' and then replace f{z) by — ln(z — z^) in the left hand 
side of ([55|) . Comparison of the net result for these two expressions yields (— ln(z — z^/), — ln(z — Zk))D 7^ 
(— ln(z — Zfe), — ln(z — Zfc'))|). The difficulty here is that In — Zfc] is not in the class of admissible analytic 
functions as specified by Not only does lim|2,|^o2 ln(z — z^) 7^ 0, but this limit is not even bounded. 
Likewise, of course, / ^ [x — XkY + (y ~ VkY] is not in the class of admissible functions for M^. The 
solution is simple. Consider the complex case first. Since the sum of two analytic functions is analytic, it is 
obvious that 

ik{z) In - In = In {^—^\ (56) 



Zk z- Zu \Z~Zk 
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if,^k)^^T;nPl) ■ (59) 



satisfies the conditions It is thus suitable as a candidate basis function for logarithmic fits. Here zj^ is 
called the paired point of Zk and it is an arbritrary complex number, except for the restrictions that |z^| < 1 
and z'j. 7^ Zk for all k. With these modifications ([55)) can then be employed to obtain 

{f,^k)^^lr{pi)~lriip'kr) , (57) 

where (p^)* is the conjugate involution of the paired point of Zk and is thus specified by taking the conjugate 
ofp^ := Equation ([57| is self-consistent and it can be used to directly verify that {^k, £,k')*D — {£.k' , ^fe)_D 

as expected. It is clear that ^k can be used to perform various types of fits in the same way that the basis 
functions l/(z — Zk) were used; however, the choice of paired points clearly affects the resulting quality of the 
fit. For example, when the values of Zk are determined by a NLLSQ process it would be advantageous to also 
determine the paired point locations as part of the same NLLSQ process, although this clearly complicates 
the implementation. Even for a LLSQ fit, the introduction of paired points is often inconvenient since it 
effectively doubles the number of point logarithmic sources and the amount of point evaluation information 
that is required for /. Thus generally it is convenient to consider another set of candidate basis functions 
given by 

Mz) ■■= In — In i = In . ^ . , (58) 

Z-Zk Z (Z-Zk) 

where all the paired points have effectively been combined and situated at the origin. When z^ := 0, the 
paired conjugate involution points are at infinity and the second term on the right hand side of ([57)1 vanishes 
for all k, by convention, and thus 

1 
2 

Clearly one can handle logarithmic source fits to analytic functions just as before by using the basis functions 
^k or ipk'i however, the branch cuts associated with these logarithmic forms raise certain issues — especially 
when f{z) does not have them. Notice that these basis functions are not translation invariant whereas the 
basis functions £,k{z) specified by ([56)1 are. 

When logarithmic fits to a harmonic function W{x, y) in are considered these branch-cut concerns 
clearly no longer apply, but the requirement that the basis functions tail off quickly enough to zero at infinity 
still remains. Keeping this concern in view, consider the following set of basis functions 

= In ^ ^ ^ In ^ = In ^^^L , (60) 

\X-Xk\ \X\ \X-Xk\ 

where X := (x, yY' and Xk '■— {xk, yk)'^ satisfy \X\ > 1 and < \Xk\ < 1. Let U = U{X) denote the linear 
superposition of these logarithmic basis terms which will be used to approximate W, then 

{/(X) = ^mfe*fe(X), (61) 
fe=i 

where ruk G M. 

On comparing ([60)) and ([58)) it is obvious that the standard completion of '^k{X) is ipk{z) and thus that 
the standard completion of U is the analytic function fu {z) where 

fu{z) -.^^^iki^kiz) (62) 
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and, in general, iik G C. In common practice harmonic problems are often mapped to C and back using 
standard harmonic completion to take advantage of various techniques available in the complex domain 
(such as conformal mapping). In the present context, when this "round-trip mapping" is considered for 
actual applications several interpretational issues immediately arise in addition to those associated with 
branch cuts. First there is the minor one involving the meaning of complex values of fik- Leaving this issue 
aside, another more fundamental one arises from the act of standard completion of W itself, since there is 
insufficient given information to insure that this completion is unique. In particular, a DIDACKS fit utilizes 
only information about W at specified points in R'^ and this by itself is obviously insufficient to uniquely 
specify the form of W over the whole region of interest so that various harmonic completions are theoretically 
possible. Specifically, at the specified data points diverse analytic functions can have identical real parts 
that match W. While a minimum principle clearly underlies this standard completion process for DACKS 
fits, it is unclear whether there may be some sort of unexpected interaction with this minimum principle 
involving the possibility of complex values of and branch cuts. It is thus at least conceptually desirable 
to quantitatively analyze this problem. This can be best accomplished by proving that the results of a fit 
performed in the complex setting must match the same fit performed in the real setting, but to do this 
requires that a unique corresponding real setting theory be available. Moreover, if this real setting theory 
exists it far easier in both a practical and theoretical sense to simply use it and side-step this round-trip 
mapping altogether. Fortunately, developing a corresponding theory based solely on the real Dirichlet form 
D[/, g] is fairly easy. The basic strategy will thus be to map the underlying inner products themselves from 
C to R^ (which only needs to be done once) rather than map the functions from R^ to C and back for 
every application. (This is clearly the most desirable way to proceed in many other situations encountered 
in complex variable/R^ potential theory as well.) 

Let G{X) and H{X) be two functions as specified from the context of ([55)1 and let h{z) and g{z) be their 
corresponding standard completions. Substituting 

dg_ _dG_ ^.dG_ (63) 
dz d X dy dz dx dy 

D[G, i/] -Re{Dj5, /i]} (64) 

{G,H)^ = Re{{g,h)^} , (65) 

{G,H)^:=^Y}[G,H] (66) 

defines the energy inner product. Here, aside from a units dependent constant, the energy associated with 
the field E = -VW is given hy \\W \\% := {W, W)e ■ In order to perform various types of DIDACKS 
fits, all that is left to do is to determine closed form expressions for (G, Kd)e, where Kn — KoiX, Xk) is 
any suitable DIDACK. Using ([55)) . it is clear that most of the work has already been done in the form of 
(Pl]) . ([55)1 and ([55]) . For example, since ipk{z) is the standard completion of 'l'fe(X), ([55]) immediately yields 

(G, vl/,)^ = iRe{(g(pD)*}- (67) 

Further, recalling that G{X) — Re{g{z)} allows for the reinterpretation of the right hand side of (|57|) . but 
the occurrence of here is perhaps disconcerting at first. Notice however that argjzfe} = arg{p^} = Ok, so 
that in the polar representation Re{(g(p^)*} ~ G(l/r/j, 6k)- Introducing Pk :— Xk/\Xk\'^ allows (|57)) to be 
reexpressed as 

iG,\ni\X\/\X-Xk\))E^lG{Pk) . (68) 



into ([39]) shows that 

and thus 

where 
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Similar higher order pole expressions can be obtained directly from (j68p by taking partials on both sides 
with respect to xu or and recalling, for example, that 



d 1 



d 1 



(69) 



dxk \X-Xk\ 



dx \X-Xk\ 



Basing fits on the expansion (j6ip and it generalizations is thus straightforward since partials of '■= 
\\W — U\W with respect to rrifc (or other source terms) are be easily determined and, as before, when the 
results are set to zero closed-form linear equation sets result that can be easily solved. Obviously ([55)1 implies 
that fits based on specific values of VW at selected points are dipole fits; moreover, similar correspondences 
between higher order poles (quadrupoles, etc.) and their associated point measurement values hold. 

Various other possibilities clearly exist besides those outlined in the preceeding sections. For example, in 
the complex setting, based on the fact that l/{z — Zk) is a shared replicating kernel for the standard integral 
norm and D norm, an alternative DACKS norm can be formed by taking a linear combination of these two 
norms: (/, g)D + A(/, g)^ , where A is real and positive. Another possibility arises from the fact that a 
single norm may have two (or more) different replicating kernels: thus in the complex setting an alternative 
DIDACK arises from the linear superposition of {z — Zk)~^ and ln[z/(z — Zk)] [c.f., (j44|) . ([58|l and ([59|l ]. 
Although this and other such possibilities are interesting in themselves, they will not been addressed in the 
sequel. 

8 Numerical Tests 

This section presents numerical test results from various SIDACKS and DIDACKS implementations. The 
emphasis here is on the complex setting and, in fact, real results will not be considered until the occurrence 
of (ffTj) towards the end of the section. 

For each implementation there are three primary aspects that must be considered: 

1. The choice of basis functions. 

2. The choice of norm. 

3. The nature of the function to be approximated. 

The third item on this list will be considered in more detail after computational issues associated with the 
first two items are fully addressed. Thus while the choices of basis functions and norms considered here, as 
well as their merits and the various issues involved, have been previously discussed and thus it is generally 
only necessary to take note of particular usages in the sequel and then observe the quality of the resulting 
fits; there are, however, associated specific computational issues that merit further discussion. In particular, 
the numerical nature of the linear equation set to be solved is clearly depend on the choice of norm and basis 
functions. Since, in all cases, (O characterizes the resulting linear equation sets, it is useful to have some 
suitable measure of the computational stability of this equation set or of the inverse of the matrix T that 
occurs there. The condition number of T is used for this purpose and it is taken to be the ratio of the largest 
to the smallest eigenvalue of T. (Since all of the eigenvalues encountered here are positive real numbers, the 
absolute value need not be taken here.) 

It is worth noting in the sequel that while the results are good, condition numbers for particular test cases 
are often surprisingly large. As discussed below all three of the factors mentioned above have a bearing on 
this phenomenon, but the effects due to basis functions clearly dominate since with a choice of orthonormal 
basis functions the computationally trivial situation T = I results (where I is the identity matrix and the 
condition mrniber is thus 1). This is one of three interconnected aspects to be considered in conjunction 



31 



with the basis functions dealt with here: their spatial localizability, their frequency spread and their non- 
orthogonality. By way of concrete example, consider the basis functions \/{z — Zk) in conjuction with the 
SIDACKS treatment of Section [51 As \zk\ approaches 1 this basis function is clearly suitable for matching 
very localized effects that occur on the unit circle, but in more common situations non-localized effects must 
be matched. In this sense the properties of the function to be approximated matter a great deal: if only a 
few non-overlapping localized effects on the unit circle are to be matched then a very well conditioned set 
of sources with \zk\ ~ 1 can be utilized. It is primarily when there is a sizeable low-frequency content in 
the function to be approximated that large condition numbers are naturally encountered. From a frequency 
perspective, when specialized to the case r — 1, pT|) yields a Fourier- like series that has coefficients with a 
magnitude of [c.f., ([6])] — which is clearly very spread out in the frequency domain unless \zk\ is small. 

While it is a property of low frequency functions to be spread out in the spatial domain, basis functions of the 
form l/{z— Zk) have a great deal of spatial overlap for \zk\ ~ and this directly implies a high degree of non- 
orthogonality or intercorrelatedness that, in turn, is responsible for the large condition numbers encountered 
here. As discussed at the end of Section SI the norms themselves also play a role here since they fail to 
exert a direct control over the source coefficients themselves. The reasons for this are linked to underlying 
aspects of the modeling philosophy itself and are worth emphasizing. For modeling problems the criteria for 
success is fairly direct and is built explicitly into the definition of the norms themselves and it is these norms 
that are minimized. Thus, while the error in the coefficients associated with a particular fit may not be well 
controlled, this generally does not matter since these errors do not show up in the norms. Moreover, for a 
pure modeling problem it is easy to validate a candidate fit at any number of points and directly ascertain 
its suitability. Hence for modeling problems condition number considerations only matter when they directly 
effect the quality of the fit itself and this often occurs only as a second-order effect when machine round-off 
effects become significant — provided, of course, that the sources are located at appropriate points. 

Several specific computational details with regards to the solution of ([5]) and condition numbers are worth 
addressing. Obviously, to obtain reasonable results in the presence of large condition numbers some care 
with regards to computational word length and equation solving techniques are required. Toward this end all 
standard complex variable computations were performed using a 64-bit word length, which corresponds to 
about 16 significant digits on the system used (on many systems 64-bit words have one or two less significant 
digits). These complex number evaluations were then mapped to and from real linear equation sets via 
(I16p . which results in each eigenvalue being duplicated (as discussed below). These real linear equations 
themselves were solved using Householder triangulation implemented with a 128-bit word length (33-34 
significant digits). (If the system architecture makes use of the fact that part of the original word length 
is reserved for the mantissa — as the system employed does — then it is possible to obtain more than double 
the number of significant digits by doubling the bit size.) The eigenvalues and thus the condition numbers 
themselves were computed using singular value decomposition (SVD) software implemented with the same 
128-bit word length. (All the computed eigenvalues encountered were positive as expected.) 

When testing the modeling application algorithms proposed here, it is useful to consider the behavior 
that might be encountered with regards to not only various basis functions and norms, but also with regards 
to various types of functions that might be fit. First, the specified class of admissible functions — which 
is somewhat restricted — must always be kept in mind. This is not a particular issue with respect to basis 
functions since all of the basis functions considered here satisfy the needed requirements, but since / is usually 
an externally specified function, it is not necessarily a given that it can be recast in a suitable form, even 
after extensive modifications. Consider, for example, the form f{z) — 1/y/z — V r^^e^*^/^. Here (/(z) — aq) 
cannot be expressed in the form of ([9]) even though it satisfies the Cauchy-Riemann conditions in polar form 
for all 1 < l^l < cxD, since it is obvious that l/\/z is not analytic at \z\ = oo, just as ^/z is not analytic at 
z = 0. (To the contrary, suppose that a series of the form 1/y/z = X]^o '^^ exists, then squaring both 
sides of this form and identifying like powers of 1/z yields an inconsistent result.) Clearly the requirement 
that candidate functions to be fit satisfy the demands of ^ is somewhat restrictive and, in fact, many 
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Function Type 


/i 


/2 


/3 


h 


h 


/6 


Average Magnitude 


0.87 


0.97 


0.90 


1.01 


1.00 


1.01 


Maximum Magnitude 


2.23 


1.38 


1.08 


1.18 


1.09 


1.72 


II II. 


1.07 


0.99 


0.91 


1.01 


1.00 


1.13 



Table 1: Summary Data for Specified Functions on Unit Circle 



standard analytic functions defined in the exterior of a unit disk do not satisfy this criteria. For example, 
sin(2;) is obviously analytic and bounded as \z\ — > oo; however, it does not satisfy the necessary criteria since 
lim \z\^oo fiz) fails to hold; moreover, even subtraction of a constant fails to produce an acceptable 
function. Clearly when a LLSQ SIDACKS or DIDACKS fit is performed to a function / that is not in the 
class of admissible functions then large errors will (generally) result since, although the replication properties 
at the data points are enforced, the minimum norm process will be satisfied not with regards to the specified 
/ as intended, but rather with respect to some function that satisfies the same replication constraints and 
that is in the class of admissible functions. In such cases some other approximating procedure should clearly 
be considered. Nevertheless, it is not difficult to come up with suitable functional forms for testing that are 
based on standard analytic functions, which may not themselves be in the class of acceptable functions. The 
last three of the six test functions to be fit were chosen with this consideration in mind. Since the discussion 
of the three items above is complete, the full list of test functions and other implementation detail will be 
considered next. 

Fits of various types were performed to one of six specified test fmictions, denoted fj for j ~ 1, 2, 3, 4, 5, 6 
and given by: 

• fi ■■= Ci E 2-" where Ci 

• /2 C2 E where C2 (1 + «)/6 and z„ := (3/4) e^'^("-i)/2 

n—1 " 

• /3 C3 E In [z/{z - z„)] where C3 := 2(1 + i) and z„ := (3/4) e^-("-i)/2 

n=l 

• /4 := sin(l/z) 

. h 2[cos(l/z) - 1] 

• /e - 1 

The scale of each of these test functions was chosen so that the resulting approximation errors would be 
more-or-less comparable. Table [T] summarizes data that gives a feeling for the properties of these various 
test functions. It was obtained by computing the magnitude of each /j at 1,000 uniformly spaced points on 
the unit circle. 

For testing purposes single ring configurations of a specified source type were used as basis function sets 
and, as mentioned earlier, only the complex setting was considered first. (These choices were made primarily 
for convenience.) The basis function radius, is shown for various specific test case configurations in 
Table [2j The basis functions used were either Xj^z — Zk) or -0^ [c.f., (|58p ]. as also indicated in Tabled For 
both sets of basis functions their locations Zk were specified by 

Zk := Rb e*2^('=-i)/^'' where i?B e M and < i?s < 1- (70) 



33 



CclSC 

Number 


Norm. 


1 n 1 nn 

to be Fit 


Basis Function Specification 


r nn pI"! nn 




k 


A- j f-j A-KjL 111 . 


1 


u 


h 


i/[Z — Zk) 


U.D 


Q 

O 


INO 




z 


u 


h 


1/[Z-Zkj 


U.D 


Id 


iNO 


o 
o 


u 


£ 

h 


1 /f^ -, \ 
L/(Z — Zk) 


1 /O 

1/2 


Id 


INO 


A 

4 


L) 


£ 

h 


1 /f^ -, \ 
L/[z~Zk) 


V.Zi 


Id 


iNO 


c 



a 


£ 

h 


l/[Z-Zk) 


1 /o 

1/2 


Id 




D 


u 


£ 

h 


Wk 


1 /o 
1/2 


Id 




i 


u 


£ 

h 


1/(Z-Zk) 


1 /o 
1/2 


Id 


AT^ 


Q 
O 


u 


£ 

h 


1/[Z-Zk) 


(J. ( 


Id 


AT^ 


n 

y 


u 


£ 

h 


1/(Z-Zk) 


1 /o 

1/2 


Id 


AT^ 


iU 


u 


£ 

h 


1/[Z-Zk) 


(J. / 


Id 


AT^ 


il 


L) 


£ 

h 


Vk 


1 /o 
1/2 


Id 


AT^ 
iNO 


1 


L) 


£ 


Ipk 


U. ( 


1 A 
10 


iNO 




n 
u 


£ 

h 


^k 


1/2 


16 


No 


14 


D 
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17 
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D 


h 
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17 
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Table 2: Case Number Specifications 
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Here Nk, the number of basis functions used, is also indicated in Table [H as is whether the complex 
energy DIDACKS norm (D) or the SIDACKS norm (ct) was used. Finally for each case the function to be 
approximated is also indicated in Table [5] so that it contains all of the relevant data for the test configurations 
and associated case numbers, but it does not contain the results. The results for each of the cases specified 
in Table [2] is presented in Table [3l The results of each fit were evaluated at 1,000 uniformly spaced points 
on the unit circle {Re = 1) and on a circle of twice that radius {Re = 2). 

A comparison of the case descriptions for cases 1 & 2 given in Table [2] with the results specified by Table [3] 
indicates that when the number of fitting points {Nk) is increased the accuracy improves, as expected. (If 
the added points are at irregular and inappropriate locations, this well may not be the case — especially for 
large condition numbers.) Comparison of the results for cases 2 through 4 indicate that as Rb becomes 
smaller, the error decreases — but there are limits. (Below Rb = 0.27 the standard deviation may be smaller, 
but not consistently so — indicating the onset of unwanted machine round-off effects.) Comparison of cases 3 
& 5 indicates that the a and D norms give roughly the same accuracy at Re = 1- Comparison of cases 3, 
7, 9, 15, 17 and 19 shows that the approach works well for various types of more or less smooth functions. 
A comparison of cases 6, 11, 13, 16, 18 and 20 shows that the same conclusion holds for logarithmic basis 
functions as well as for simple poles. As expected from the above discussion, comparison of cases 7 & 8 
and/or 9 & 10 shows that smaller values of Rb may not be appropriate for certain functions and, in fact, 
that the choice of source depth {Rb) and spacing (i.e., Nk) may be highly dependent on the nature of the 
function to be fit. Cases 11 through 14 show that this conclusion is not particularly dependent on the choice 
of basis function type. In this connection, notice that in some sense /2 and fs are more irregular than /i 
since the series representations of /2 and /s in the form ([9]) do not terminate. Functions f^, and /g are 
also very well behaved and the effects of basis function choice for each of these function types can be had 
through a study of cases 14 through 20. Finally the inclusion of a single additional pole source at the origin 
(i.e., 1/z) improves accuracy, but raises the condition number, as is apparent from comparison of cases 21 
through 23 with cases 5, 3 and 6 respectively. Computationally, for the a norm these 1/z terms are easy to 
include since it is only necessary to add an appropriate row and column to the matrix ([5]). For example, a 
direct substitution of 1/z and the series representation lU into (fT2| yields the needed additional rows and 
columns to fill out the T matrix: {1/z, {z — Zk)^^)a ~ {{z — Zk)^^, l/z)a = {1/z, l/z)„ — 1. Also when a 
row is added to A the additional element is A]^f,^\ = a\. The various required terms for the D norm are 
obtained in a like fashion. 

Obviously the condition numbers given in Table [3] are somewhat large. Since the matrix T occuring 
in (|15p or (|47p is Hermitian, the associated eigenvalues are real. Further, since the norm of any function 
here is positive semidefinite and exists by assumption all of the T matrices encountered in testing were 
positive, it follows that |T| = \k > holds, where the Afc are the Nk eigenvalues of T. As previously 
mentioned, each of the eigenvalues Afe is repeated twice for the real symmetric matrix occuring on the left 
hand side of and thus, when T is represented in this form, |T| can be computed by taking the square 
root of the product of all of its eigenvalues. Since a SVD analysis of the symmetric real valued form of T was 
performed in all cases, the results of this determinant computation were readily available. For SIDACKS 
cases this square root computation can be directly compared to the complex number evaluation of the right 
hand side of (fT9|) . This comparison was performed for all the relevant cases in Tables [2] and [3] and the results 
matched to within the expected accuracy. For example. Case 5 yielded .1044048746849805x10"^^ for the 
square root eigenvalue computation versus .1044048714876914949x10"^^- .52161 x 10"*^* z from (fT9)) . 

Finally, limited testing based on ([60| and ((68)) was performed in the real setting. Similar results to 
those found in Table [3] were obtained. Results from a single case here will have to suffice. The function to 
be approximated was defined as 




4 



(71) 
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Oase 
Number 


Condition 
Number 


Fit Errors at i?_E = 1 


Fit Errors at = 2 


Standard 
Deviation 


Maximum 
Magnitude 


Standard 
Deviation 


Maximum 
Magnitude 


1 


.160x10^ 


.180x10-1 


.377x10-1 


.916x10-3 


.116x10-2 


2 


.283x106 


.302x10-3 


.643x10-3 


.488x10-6 


.605x10-6 


3 


.671x10^ 


.163x10-4 


.348x10-4 


.143x10-® 


.189x10-® 


4 


.701x1016 


.123x10-® 


.223x10-® 


.105x10-12 


.118x10-12 


5 


.107x10^" 


.163x10-'' 


.349x10-4 


.125x10-9 


.257x10-9 


6 


.172x10" 


.214x10-5 


.469x10-5 


.881x10-11 


.225x10-1° 


7 


.671x10® 


.995x10-2 


.138x10-1 


.123x10-5 


.130x10-5 


8 


.277x10"^ 


.667x10-2 


.889x10-2 


.179x10-3 


.179x10-3 


9 


.671x10^ 


.185x10-2 


.243x10-2 


.866x10-® 


.103x10-^ 


10 


.277x10^ 


.118x10-2 


.122x10-2 


.123x10-5 


.123x10-5 


11 


.172x10" 


.997x10-2 


.138x10-1 


.102x10-6 


.146x10-6 


12 


.710x10^ 


.971x10-2 


.132x10-1 


.154x10-4 


.157x10-4 


13 


.172x10" 


.186x10-2 


.244x10-2 


.242x10-® 


.349x10-® 


14 


.710x10^ 


.124x10-2 


.163x10-2 


.249x10-6 


.253x10-6 


15 


.671x10^ 


.155x10-4 


.179x10-4 


.198x10-® 


.213x10-® 


16 


.172x10" 


.984x10-6 


.133x10-5 


.975x10-11 


.151x10-1° 


17 


.671x10** 


.153x10-4 


.166x10-4 


.527x10-9 


.589x10-9 


18 


.172x10" 


.172x10-'^ 


.196x10-5 


.916x10-11 


.134x10-1° 


19 


.671x10® 


.173x10-4 


.262x10-4 


.200x10-® 


.242x10-® 


20 


.172x10" 


.131x10-5 


.231x10-5 


.108x10-1° 


.218x10-1° 


21 


.776x10" 


.122x10-4 


.241x10-4 


.452x10-1° 


.955x10-1° 


22 


.456x10^" 


.122x10-4 


.241x10-4 


.276x10-9 


.417x10-9 


23 


.456x101'' 


.204x10-5 


.406x10-5 


.554x10-11 


.131x10-1° 



Table 3: Fit Results for Cases Specified in Table [2 
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where X'^^ := ^(cosfJ,^ , sinJ^J-^ and where (5'j :— {n — 1)tt/2. When evaluated at r = 1 this specified function 
definition yielded Max{ \F{X)\} = 1.785 and Avg{|F(jf)|} = 1.130. A real linear combination of 16 
basis functions specified by (j60p were used as an approximating form, where Xk ■= (2/5) (cosAfc, sinAfc)-^ 
and Afc {n — 1)tt/8. The condition number for the resulting T matrix was .6939 x 10^^. Evaluated at 
Re = 1, the resulting magnitude of the fit errors had a standard deviation of .4124 x 10"'' and a maximum 
of .6053 X 10~^. The corresponding errors at Re = 2 were a standard deviation of .4342 x 10—7 and a 
maximum of = .759784 x 10^^. Observe that since the fitting parameters are real there are half as many 
parameters as there are for the typical fits given in Table [3l Although the condition number was large here, 
the smallest eigenvalue was an outlier and excluding it in the inversion process (which can be easily done with 
Householder or SVD software) produced a fit error with almost identical errors that had a condition number 
of .1490 X lO''. While this is obviously one way to circumvent high condition numbers, it is also interesting 
that due to the smoothness of F and the large word length (128-bit) good fits here can be obtained even for 
small values of Rb (even below .01). (Using Rb = -01 in the above basis function specifications and fitting 
to the same F yields a fit with a standard deviation of .3834 x 10"'' and a maximum of = .5469 x 10~^ at 
Re — 1. The condition number was .7517 x 10"^^.) 



9 DACKS Theory for Disk Interiors 

Finally SIDACKS and DIDACKS theory so far has focused exclusively on the exterior of a unit disk, but 
these DACKS approaches can easily be recast in terms of unit disk interiors. For convenience, it is useful 
to reuse the same symbols as much as possible, so it will be assumed in the sequel that oo > |zfc| > 1 > |z| 
specifies the basic parameters of interest for the interior case. With this understanding, the inner product 
specified by (jl2p can be utilized as is with an appended subscript a' to replace the a subscript — where, as 
usual, the prime will be used to connote a complimentary region. (While the compliment of the unit circle 
is the unit circle, the analytic region of interest is different for a and a' .) Further, it is useful to specify the 
admissible analytic functions for this interior geometry by upper case letters rather than lower case letters 
in order to keep in mind that they have a different series representation. Thus consider a given function 
F(z) that is to be approximated by a linear combination of simple poles located in some proper subset of 
the disk's exterior and which have the form l/{zk — z). (The reverse sign convention is taken here so that 
l/(zfe — z) will be positive when Zk and z are located along the positive real ai-axis.) In the present context 
([9|) and (fTTj) are replaced by 

oo oo ^ 

F{z) = Y.bnZ- and =^^. (72) 

n— n— 

Clearly, if closed-form inner products of the form {F, (zj. — z)~^)a-' can be found the rest of the interior 
SIDACKS development is straightforward. Towards that end, replacing / and g by F and l/(zfe — z) in 
and replicating the steps that led up to (fTi]) yields: 

(F, (zfe - z)-i)^, = pk {F{pl)r , (73) 

where pk ■= 1/zk as before. Recall that a discussion of the merits of the SIDACKS procedure just outlined 
and Szego kernel based interpolation schemes was given at the end of Section \5\ 

DIDACKS theory can be also be developed for the interior of a disk in much the same way. Using the 
previous notational conventions in an obvious way, the energy inner product for the interior of a unit disk 
can be defined as 

r=0 9=0 
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Now, however, F and G are assumed to have a power series expansion of the form 



f^(z) = ^6„z", (75) 



n=l 

where the constant term (60) has been omitted from the analogous series given in ([72|) . This constant term 
was also omitted from ^ for reasons discussed in Section [H This means that 

Uz) := - - = , ' , = E ^ (76) 

will be used as basis function here instead of l/{zk — z). Replacing G in (|74p with this series expansion and 
replicating the steps that led to (|44p yields the following result 

{F,{z-Zk)-')^,^^{FM)r ■ (77) 

Interior complex SIDACKS and DIDACKS theory can obviously be easily built upon the above results in 
complete analogy with the exterior treatments in Sections O and [HI thus it is not necessary to include the 
details here. The development also proceeds in a completely analogous fashion to that of Section [71 
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Appendix A 



Discussion of Educational Aspects 

As noted in the introduction, there are various educational possibihties that are inherent in the general 
DIDACKS approach. Since these education issues will not be explored systematically in any other venue by 
the author, several relevant points are summarized here. After these points are made, this appendix ends 
with a brief discussion of self-educational aspects that, as indicated in the introduction, are considered to be 
of equal or even greater significance than these educational ones. 

Before considering these issues it is useful to take a step back and briefly look at the bigger picture. While 
much of the article has been written at an accessible level, this still does not mean that it was intended as an 
introductory one and, as such, it does not really stand alone in this regard (i.e., neophytes probably should 
have the guidance of a more seasoned intermediary before approaching much of the material here). In order 
to simplify matters in general, although the author might have wished it differently, with regards to the 
overall level of material there is a curious attribute of DIDACKS theory itself that it has an certain innate 
multi-scale accessibility, which is to say that any article, even a relatively brief one (which this one is not) 
can, in some sense or other, cover a very large span of readership — from material that is at the introductory 
level to material dealing with open research problems; moreover, at both ends of the spectrum the problems 
and challenges posed seem to be unique and interesting. For example, at the introductory or educational 
level, it is possible to contemplate combining the material of many core courses into one DIDACKS centered 
course, so that material that might not otherwise be covered in a core curriculum can be covered. Related 
educational items are: 

1. DIDACKS theory can serve as a means of introducing upper level undergraduate or graduate students to 
approximation theory, harmonic analysis, functional analysis, integral equation theory, computational 
linear algebra, kernel theory, Bergman kernel theory, RKHS theory and other advanced topics. Students 
who have some previous exposure to complex variables will be familiar with the functional form l/{z — 
Zk) and it is not too big a conceptual jump to consider it as a kernel within the wider context of 
DIDACKS theory. In particular, RKHS theory should be considered for inclusion in introductory 
functional analysis courses, but it usually is not, in part, due to conceptual hurdles, and DIDACKS 
theory may be a way of lowering these hurdles. One way to overcome these conceptual hurdles is to 
have students take a hands-on approach by having them develop thir own software implementations. In 
this regard, the reproducing property of RKHS fits allows the students to validate their own software. 
The basic DIDACKS approach is even more ideally suited to this style of supplemented self-instruction 
than RKHS theory is due to the overall accessibility of the material. 

2. Besides trying to motivate others to consider injecting RKHS theory, where appropriate, into classes, 
one goal here that is even more central is to motivate lecturers to present discussions of Dirichlet's 
principle and of the associated Dirichlet form (or integral). Others have suggested that Dirichlet's 
principle should be a main part of various introductory mathematics courses due to its central role 
in the history of various branches of analysis. One existing significant realization of this plan is [7], 
which covers real analysis via variational principles. The material covered here is another way to inject 
Dirichlet's integral and Dirichlet's principle into the coverage of introductory material. 

3. Aside from the natural way that other areas of mathematics can be integrated into a consideration 
of the DIDACKS approach, there are also a very wide range of tangentially related scientific and 
engineering disciplines that it can be applied to or adapted to in a non-trivial way. By this means 
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applied mathematicians, or other students in question, can be given a very broad exposure to diverse 
areas. 

4. There are undoubtedly pedagogical advantages to taking a single-topic or focused approach to the 
coverage of basic material. One advantage is the possibly of an increased level of interest. A second 
is that students often mentally retain mathematical tools by distinct separate categories, which are 
linked to only specific subjects, and thus miss common mathematical themes; hence, studying topics 
from the perspective of a different coherent mathematical framework could help prevent this. Along 
these same lines, with a coherent theme it is reasonable to expect that students will do a better job of 
integrating mathematical material with applications and thus retain it better. From a research point 
of view, these same attributes clearly confer cross-fertilization possibilities. 

5. In connection with point 1 above, it is also possible to use the DIDACKS approach as a framework 
for a self-instruction courses in scientific or applied programming courses where numerical consistency 
checks are extensively utilized. For example, in the complex setting, for some specified / by taking 
finite numerical differences along different directions and at different points, the fact that df/dz is 
both computed correctly and that / is analytic can be easily directly spot-checked by the student 
programmer. While the primary implementation tool for overall consistency checks is the replication 
property associated with the DIDACKS interpolation fits themselves, it should be noted that symmetric 
configurations of interpolation point positions in both the complex and real plane may well produce a 
consistent collocation check result even though implementation errors are present due to the existence 
of various underlying reflection and symmetry principles. Thus a software consistency check should 
entail using not only irregular interpolation point locations, but also an / that varies in value from one 
interpolation point to the next. 

6. Section O and Section [5] are largely self contained. Since the material in Section [5] up through (pi]) 
does not require a significant functional analysis background and is accessible at a very low level, when 
placed in a suitable context it can be used as a stepping stone to more advanced topics; moreover, it 
is possible to design a course around the content of Sections [5] and [6] that introduces a variety of pure 
and applied concepts. Furthermore, after closer examination of the content of this article, it should 
be apparent that much of the material often present in an introductory course in complex analysis 
has been covered and that those topics that have not been covered can be introduced with a little 
ingenuity. For example, one topic not covered is conformal mapping, but it is obvious that this topic 
can be naturally integrated by introducing pre-mappings from other regions to the unit disk (and then 
return post-mappings). (Other classes of functions besides analytic functions can even be considered in 
this connection provided that the class of mappings entertained here is sufficiently broad.) One other 
common topic not directly dealt with in the body of the paper is residue theory — see the next item. 

7. Those students who are already knowledgeable in the theory of residues and poles, as well as some 
others, may slip into a momentary lapse of confusion and think that since residue theory can be used 
to evaluate isolated poles at known positions by performing a simple closed line integral, this process 
can be used here too and thus that the formalism developed herein is, in some sense, superfluous. To 
see the issue here, consider the case where the functions in question, /(z), are analytic in the exterior 
region and all their poles are located in the interior region at unknown locations. The point, of course, 
is that it is assumed that noting is known even in principle about the values of f{z) for |z| < 1 and all 
that is known in this interior region is the location of the poles that are to be used in the approximating 
form; hence, any part of a line integral of / that is inside the unit disk is unknown. (As discussed 
below, / may well arise from a continuous source distribution rather than from a denumerable set of 
isolated poles.) By directly trying to use the theory of residues under these restrictions the student can 
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convince him or herself that residue theory cannot otherwise supply a sufficient amount of information 
to uniquely determine the pole strengths. 

From another perspective, consider the question as to whether some combination of simple poles in the 
interior of a disk exists that can produce a null result in the exterior region: /(z) = for |z| > 1. This 
question cannot be directly resolved by using line integrals in the exterior region, but it is answered in 
the main body of the text by other techniques. 

8. While it may not at first be obvious, with a sufficient application of ingenuity numerous and varied 
interesting problem sets can be generated to accompany lecture notes based on the DIDACKS approach. 
One example here will have to suffice: Find a polynomial expression in terms of u{9) = e'^ for the 
extremals of [given by (fTU]) ] on the unit circle. (Notice that for Nk < 4 closed-form expressions can 
be written out for the pole strengths.) Solution: (a) Take d /dO of (p*(p and set the result to zero, (b) 
Substitute 1/u for u*. (c) Multiply through by u and then by the common denominator. 

9. A thorough discussion of the interplay of power series and analytic functions is rightly considered to be 
a core part of standard complex variable courses. In conjunction with the precise minimal requirements 
here for series or of the corresponding smoothness requirements for harmonic functions, interested 
readers are encouraged to consult standard complex variable or harmonic analysis texts [2] and rethink 
the associated DIDACKS requirements for themselves. With regards to the discussion of these partic- 
ular topics elsewhere in this article, a certain amount of looseness was present so that instructors who 
choose to present the associated material can discuss smoothness and power series related topics to 
the level of detail and rigor that they wish. One consideration here that is important for DIDACKS 
theory is that continuous source distributions can be entertained, which brings in concepts not usually 
addressed at the introductory complex variable level since it involves issues that are not raised when 
only denumerable sets of simple poles are considered. 

If lecturers do bring up the topics of power series in this context, certain obvious examples should also 
probably be included in the discussion, such as the harmonic problem for the exterior of a unit 
disk (fi) where / = 1 on the upper half of the unit circle and f — along the lower half of the unit 
circle. (The point here is that although the Poisson integral gives a harmonic solution /, which is 
bounded, it is not altogether obvious that D[/, /] itself is bounded.) To further explore this example 
the series solution obtained by standard boundary value techniques can be used to explicitly evaluate 

D[/, /, 1, n]. 

Autodidactic Aspects of Notation and Nomenclature 

The ability to engage in self-education has long been recognized as a fundamental trait of successful re- 
searchers. This is partially due to the fact that new ideas often come from exposure to other concepts 
from fields outside of ones narrow specialty. It was thus with some hope of alluding to certain autodidactic 
possibilities and also of certain cross-fertilization possibilities that the acronym DIDACKS was chosen. 

Unless a researcher who wants to broaden his or her horizons is lucky enough to have a colleague to 
interact with who is familiar with some outside area of interest, there are often significant natural barriers, 
in and of themselves, that must overcome in order to master this other area of expertise. By employing an 
overall accessible level of presentation some of the inherent natural boundaries erected by overspecialization 
can be overcome. Whether the present article is in any sense successful or not on this score, this factor 
is at least brought up here for general consideration by others. Finally, along similar lines, a much more 
than usual amount of attention has been given here to notation in the hopes that it might foster more 
consideration and discussion in the open literature of this particular topic by others. 
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Appendix B 



Dirichlet Integrals for Analytic or Antianalytic Functions 



In this appendix complex valued functions are consider to be either analytic or antianalytic. Such 
functions are a subset of those that depend on both z and z. The assertion that / is analytic means that 
df /dz = or that df/dz = df /dz. Likewise, the assertion that / is antianalytic means that df /dz = 0. 
After a little thought it might seem obvious that the correct generalization of the Dirichlet integral 



dz 



dg 



dA 



(B-1) 



in the case where either df/dz = or df/dz = is the expression 
©c[/, g] := 



dl 
d z 



d z 



dz 



d_g 

d z 



dA 



where a different type face has been used for the new form for the Dirichlet integral. Let 

2c[f, g] ■■= 



dl 

dz 



dg 
d z 



dj_ 

dz 



dg_ 

d z 



(B-2) 



(B-3) 



denote the integrand here. The question addressed in this appendix is how to justify this form of 3c [/i .9] 
staring from expressions that are obviously related to Dirichlet forms. 
Consider the expression 



dj_ 

dx 



d_g_ 

dx 



djV fdg 
dyj \dy 



For g = f = u + iv, (for general C°° u and v) it is easy to show that 

rr j^i fdu\^ f dv\^ fdu\^ fdv 



or that 



Jc[/, /] dA -D[u, u]+B[v, v]. 



(B-4) 



(B-5) 



(B-6) 



hence 3^c[/, 5]/2 is a good candidate to examine as the correct integrand for the more general form of the 
complex Dirichlet integral. Since dz/dx — dz/dx = l and d z/dy — —dz/dy — i, the chain rule gives 



df df df df .df .df 

t;— = t; ^ 7r= ^i^d — = I- I — 

dx d z d z dy d z d z 

A straightforward evaluation thus gives 

J,[f, g]/2 = Uf. g\ 
as expected, which justifies definition (|B-2|) from a different perspective. 



(B-7) 



(B-8) 
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